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Abstract 

Details of the analysis of the r/vr" system studied in the reaction 7r~p -^ r]TT~p 
at 18 GeV/c are given. Separate analyses for the 27 and Tr+vr^vr'^ decay 
modes of the r/ are presented. An amplitude analysis of the data indi- 
cates the presence of interference between the a^(1320) and a J = 1 ^ 
wave between 1.2 and 1.6 GeV/c^. The phase difference between these 
waves shows phase motion not attributable solely to the a^(1320). The 
data can be fitted by interference between the a^(1320) and an exotic 
1-+ resonance with M = (1370 ±16 +^°) MeV/c^ and T = (385 ±40 
_IQ^) MeV/c^. Our results are compared with those of other experiments. 



Typeset using REVT^ 



I. INTRODUCTION 

In a previous publication [|1|, evidence was presented for an exotic meson produced in 
the reaction 

'K~p — >■ ri'K~p (1) 

at 18 GeV/c, with the decay mode rj — > 77. The purpose of this paper is to provide details of 
that analysis, to discuss additional analyses of those data, and to give a detailed comparison 
of our results with those of other experiments. We also compare those results with data from 
our experiment on Reaction (|I|) but with the r] -^ tt+tt^tt^ decay. 

The Tj-K system is particularly interesting in searching for exotic (or non-gg) mesons 
because the system has spin (J), parity (P), and charge-conjugation (C) in the sequence 
J^*" = 0^^, 1^^, 2++, 3 ^ ... for £ = 0, 1, 2, 3, ... . (Here d is the orbital angular momentum of 
the rjTT system.) Hence a resonance with an rjir decay mode with odd l is manifestly exotic. p~| 
Having isospin 1=1, such a resonance could not be a glueball {2g,3g, . . .), but it could be a 
hybrid (qqg) or a multiquark (qqqq) state. 

A. Models 



Properties of hybrids and multiquark mesons have been discussed in the framework of 
various models p-p^. Calculations based upon the MIT bag model predict [§|-§[ that an 
1 = 1, 1 ^ hybrid {qqg) will have a mass near 1.4 GeV/c^. On the other hand, the flux-tube 
model 0,§] predicts the mass of the lowest-lying hybrid state to be around 1.8 GeV/c^. 
Characteristics of bag-model 5- wave multiquark states (which would have J^ = 0^, 1^ or 
2+) have been discussed but those for a 1^ state have not. QCD sum-rule predictions ]10 



'^ A qq meson with orbital angular momentum L and total spin S must have P — (—1)^^^ and the 
neutral member of its isospin multiplet must have C ~ (— 1)^+ . A resonance with quantum number 
in the sequence J^*^ = 0^^,0^ ,1 ^,2^ ,3 ^ ■■■ does not satisfy these conditions and must be 
exotic. 



vary widely between 1.0 GeV/c^ and 2.5 GeV/c^. Recently, an analysis of the multiquark 
hybrids has been carried out, based on the diquark cluster model |]TT|; this model predicts 
a lowest-lying isovector J-^*^ = 1 ^ state at 1.39 GeV/c^ with a very narrow width (~ 8 
MeV/c^). Finally, recent lattice-gauge calculations yield mass estimates for a 1""*" hybrid in 



the range from 1.7 to 2.1 GeV/c [|12 



B. Previous Experiments 

Several experiments, prior to the publication of the E852 results [0], had studied the 
rjiT final state, and observed an enhancement in the P wave around 1400 MeV/c^ |ir3|-|T6 
However, they reached conflicting conclusions. 



The 1988 GAMS experiment |T3[ at CERN {tt~p at 100 GeV/c) claimed to find a nar- 
row enhancement in the unnatural parity exchange Pq wave, but found the natural parity 
exchange P+ wave to be "structureless" . The method of analysis and the conclusions were 
seriously disputed by some of the same authors later [|l^] . 



The 1993 VES experiment |]T4[ at Serpukhov (n^N at 37 GeV/c on a beryllium target) 
found enhancement in the natural parity exchange P+ wave and concluded that "the P+ wave 
is small but statistically significant and contains a broad bump" . They made no attempt to 
identify the "bump" with a resonance. 

The 1993 KEK experiment [|1^ (7r~p at 6.3 GeV/c) claimed that "A clear enhancement 
of the P+ wave was observed around 1.3 GeV/c^" but noted that "The phase of the P+ wave 
relative to the D^ wave shows no distinct variation with mass in the analysis region" . They 
therefore made no attempt to offer a resonance hypothesis. 

The 1994 Crystal Barrel experiment [ITH] on pp annihilation at rest concluded that their 



"TT^vr'^r] data may at most accomodate a small amount of featureless ttt] P-wave" 

The first claim for a 1~"^ exotic resonance in the rjir' channel was made by our experiment 
at BNL |jl|. In this paper we present details of the measurements and method of analysis 
used in that earlier letter publication. We note that since the publication of our letter, an 



independent confirmation of our results has come from a new measurement by the Crystal 
Barrel collaboration ||18[|. 



II. EXPERIMENTAL DETAILS 



A. E852 Apparatus 



Our data sample was collected in the first data run of E852 at the Alternating Gradi- 
ent Synchrotron (AGS) at Brookhaven National Laboratory (BNL) with the Multi-Particle 
Spectrometer (MPS) [|l^ augmented by additional detectors. A diagram of the experimen- 
tal apparatus is shown in Fig. |l]. A Cerenkov tagged vr" beam of momentum 18 GeV/c 
was incident on a one- foot long liquid hydrogen target at the center of the MPS magnet. 
The target was surrounded by a four-layer cylindrical drift chamber (TCYL) ||2^ used to 
trigger on the proton recoil of Reaction (|l]), and a 198-element cylindrical thallium- doped 
cesium iodide array (Csl) |^ capable of rejecting events with wide-angle photons. The 
downstream half of the magnet was equipped with six drift chamber modules (DCl-6) | p2|| , 
each consisting of seven planes, used for charged-particle tracking. Interspersed among these 
were: three proportional wire chambers (TPXl-3) to allow triggering on the multiplicity 
of forward tracks; a window-frame lead scintillator photon veto counter (DEA) to ensure 
photon hermeticity; a scintillation counter (CPVB) to veto forward charged tracks for neu- 
tral triggers; and a window-frame scintillation counter (CPVC) to identify charged particles 
entering the DEA. Beyond the magnet were: a newly-built drift chamber (TDX4) consisting 
of two x-planes; two scintillation counters (BV and EV) to veto non-interacting beam tracks 
and elastic scatters respectively; and a 3045-element lead glass electromagnetic calorimeter 
(LGD) pH] to detect forward photons. Further details are given elsewhere 



B. Trigger 



The trigger (see Ref. [Q for details) for Reaction (|1|) required a recoil charged particle in 
the TCYL detector and one charged particle traversing each of the first two TPX chambers. 
In addition, an electronic algorithm [^ coupling energy and position information in the 
LGD calorimeter (an "effective-mass" trigger) was utilized for the purpose of enhancing the 
fraction of ?7's relative to vr'^'s in the sample. A total of 47 million triggers of this type were 
recorded. 

C. Event Reconstruction and Selection 

Of the 47 million triggers, 47,235 events were reconstructed which were consistent with 
Reaction (|l|). These were selected by requiring: 

• topological and trigger cuts including requirements for: 

1. two photons reconstructed in the LGD; 

2. one forward track reconstructed in DCl-6; 

3. one recoil track reconstructed in TCYL; 

4. a common vertex, in a target fiducial volume, reconstructed from the charged 
tracks and the beam track; 

5. no energy deposited in the DEA detector or outside the fiducial volume of the 
LGD; 

6. the energy deposited in the Csl arrary being less than 160 MeV; 

• that the effective mass of the two photons be consistent with the rj effective mass with 
a confidence level greater than 10~^; 

• that all data come from runs which had proper functioning of the trigger processor; 



that the photons hit the LGD within a fiducial volume which excluded a 4.0-cm region 
(one block width) around the periphery of the LGD as well as a 4.0-cm wide region 
surrounding the beam hole; 

that the distance between a photon and a charged track hitting the LGD exceed 20 
cm; 

-LO < (missing mass)^ < 2.5 (GeV/c^)^; 



• a SQUAW [^ kinematic fit (requiring energy and momentum conservation) to Reac- 
tion (|1]) with a confidence level > 10%; 

• that the difference in angle A0 between the fitted proton direction and the measured 
track in TCYL be less than 8°; 

• the exclusion of events for which the vr" went through an insensitive region of TPX2; 

• the exclusion of events for which the tt~ went through a small region surrounding the 
EV/BV veto counters. (Events which had a vr" in this region were sometimes vetoed, 
probably due to Cerenkov radiation in the EV or BV light pipes). 

Shown in Table | is the effect on the data sample for each of these cuts. The last two cuts 
listed in the table were additional cuts made on the data to carry out the PWA. 

D. Experimental Acceptance 

The experimental acceptance is determined by a Monte Carlo method. Events are gen- 
erated using SAGE |2^ with peripheral production (of the form dN/d\t\ = Aexp —b\t\ with 
b = 4.0(GeV/c)~^) and with isotropic angular distributions in the Gottfried- Jackson (GJ) 
frame. (The GJ frame is a rest frame of the rjir^ system in which the z-axis is in the direc- 
tion of the beam momentum, and the y-axis is in the direction of the vector cross-product 
of the target and recoil momenta.) After adding detector simulation using GEANT |^, the 



Monte Carlo event sample is subjected to the same event-selection cuts and run through the 



same analysis as the data. A second method (called SAGEN) which did not use GEANT 
was employed as well. This method also used SAGE as the event generator, but instead of 
using GEANT, the acceptance was determined using a full detector simulation but without 
such effects as multiple scattering, pair production and secondary interactions. This second 
method allowed acceptances to be calculated much more quickly. The only differences noted 
in the amplitude analysis results (discussed below) between the two methods was in the 
number of events, since GEANT takes account of pair production and secondary interac- 
tions, whereas SAGEN does not. The average acceptances are the ratios of the generated 
events to the accepted events and are shown in Figs. ^-^ using the SAGEN method. 

The average acceptance as a function of r^vr" effective mass is shown in Fig. ^ The 
average acceptance decreases by about a factor of two over the effective mass region from 1.0 
to 2.0 GeV/c^. Average acceptances are calculated for peripheral production and isotropic 
decay as described above. 

Shown in Figs. |] and |^ is the acceptance as a function of cosO and of for various ranges 
of the riTT~ effective mass. Here 6 and (f) are the polar and azimuthal angles measured in the 
GJ frame. The polar angle is the angle between the beam direction and the rj direction in 
this frame. The inefficiency in the backward direction corresponds to slow r^'s and fast tt^'s 
in the lab. The slow ?7's lead to low energy 7's which are often produced at wide angles and 
thus miss the LGD. In some cases, fast tt^'s cause the event to be vetoed if they strike the 
EV or BV scintillation counters, leading to further inefficiency in the backward direction. 
The inefficiency in the forward direction is due to an inefficiency in detecting slow, wide- 
angle pions which can scatter in the Csl detector. The acceptance in is relatively uniform. 
There is a correlation between the energy of the 77 and at finite momentum transfer and 
this leads to the observed shape. 

Finally, shown in Fig.|^ is the average acceptance (GEANT-based) as a function of |t'| = 
\t — tminl; whcre t is the the four-momentum-transfer between the initial- and final-state 
protons and tmin is the minimum value for this quantity for a given tjit' effective mass. 
The dramatic decrease in acceptance below about |t'| = 0.08(GeV/c)^ is due to a trigger 
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requirement. In particular, since we require the presence of a recoil proton in TCYL, the 
trigger cannot be satisfied if the proton stops in the hydrogen target. 

E. Background Studies 

Shown in Fig. ^ is the 27 effective mass distribution for events in the a^(1320) mass region 
from 1.22 GeV/c^ to 1.42 GeV/c^. The data sample used for this distribution consisted of 
a subset of events satisfying the cuts listed in Table | but without SQUAW confidence level 
cuts. The central cross-hatched region in Fig. ^ shows the events which remain after the 
SQUAW-based kinematic-fitting cuts. The distribution has a ~ .03 GeV/c^. Both this 
distribution and the missing-mass-squared distribution discussed below are consistent with 
that expected from Monte Carlo studies when the energy resolution of the LCD for a photon 
of energy E is taken to be of the form a/ E = a+b/y^ with a = 0.032 and b = .096 (GeV)^/^. 
(This was the resolution function used for the LGD in the kinematic fitting.) Two methods 
have been used to study the background in our sample. Method 1 used the shaded sidebands 
of Fig. P and allows us to study the non-r^ background in the data. 

The missing-mass-squared distribution for the data sample before kinematic fitting is 
shown in Fig.^ The dashed histogram shows the events which remain after kinematic 
fitting. The distribution for good events for Reaction (|lD should peak at the square of the 
proton mass or at a value of 0.88 (GeV/c^)^ . 

A scatterplot of the missing-mass squared versus the 27 effective mass is shown in Fig. ^ 
Background studies using our Method 2 take as the background estimator a region sur- 
rounding the central signal region seen here instead of using the sidebands of Fig. |^. In this 
way we take into account background events of both the non-?7 type and of the type which 
does have an 77 present but is not exclusively Reaction |l| such as events with an extra 7r°. 
(The background is estimated using the region included within the outer elliptical area of 
Fig. H but not within the middle elliptical region. Events in this elliptical band are used to 
determine the magnitude of the background as a function of rjir effective mass.) 



Shown in Fig.|| is the effective- mass distribution of the background, estimated using 
Method 1. In this figure are shown the effective-mass distribution for each side-band region 
as well as the summed distribution for the background regions. Because the background 
regions have different thresholds, one higher than the signal region and one lower than the 
signal region, the histograms are shifted by an appropriate amount (so that their thresholds 
match that of the signal region) before summing. 

In Fig. |10| is shown the polar angular distribution of the background events from Method 
1 in the 02^(1320) effective-mass region from 1.22 GeV/c^ to 1.42 GeV/c^. Distributions are 
shown separately for the low-mass and the high-mass sidebands of the rj. This high-mass 
sideband distribution is somewhat peaked in the backwards direction with a tendency for 
the distribution to have an excess below the region cos^ < —0.5. We note that this is in the 
opposite direction from the asymmetry in the data (see below) and therefore cannot be the 
cause of the observed asymmetry. Of course the intensity of the background is quite small 



as seen in Fig.|Tl| below and therefore could not lead to significant changes in the angular 



distributions in the data in any case. 



III. GENERAL FEATURES OF THE DATA 



The 02^(1320) is the dominant feature of the t^tt" effective- mass spectrum shown in Fig, p 
The background, which is shown shaded in the figure, is estimated from Method 2 above, 
and is approximately 7% at 1.2 GeV/c^, and only 1% at 1.3 GeV/c^. 

The acceptance-corrected distribution of |t'| = \t — tmini is shown as the solid points 
in Fig.O for |t'| > 0.08(GeV/c)^. (Our acceptance is quite low below 0.08 (GeV/c)^ as 



discussed in Section [11 D| .) Since the data are dominated by a^(1320) production, we show 
as the solid curve (1) the prediction of a Regge Pole model for the differential cross section 
for the reaction n^p —* a^p at 18 GeV/c. (Note that the ordinate values are given by the 
theory and the data are normalized to the theory so we are comparing only the shape of 
the data with the theoretical prediction.) The model includes contributions from p and /2 



Regge trajectories with parameters from a fit by Sacharidis [^ and also includes a small 



/o at t' = 0.15(GeV/c)^) contribution from a uniform {t' independent) background. For 
comparison, results of another experiment (open circles) [^ which studied the reaction 
TT^p -^ a^p, a^ —* K^K^ at 22.4 GeV/c and the Sacharidis fit (curve 2) are shown. (Note 
that the values shown take into account the 02 decay branching fractions.) We conclude that 
the shape of our t' distribution is consistent with previous experiments and with natural- 
parity exchange production in Regge-pole phenomenology pO||. 



Acceptance-corrected distributions of cos 6' are shown in Fig. |T3|for various ranges of r]7r~ 
effective mass. For illustration purposes, the acceptance correction is calculated here for 
isotropic decay of the 777?" system. The acceptance correction used in the amplitude analysis 
discussed below is based upon the observed decay angular distribution. The presence of a 
significant forward-backward asymmetry in the cos^ distribution is obvious. 

The forward-backward asymmetry in cos 6 is plotted as a function of 777?" effective mass 



in Fig. |14[ Here, the asymmetry is defined as {F — B)/{F + B) where F{B) is the number 
of events in the mass bin with the t] decaying forward (backward) in the GJ frame. For this 
figure, the asymmetry was calculated for events in the region with |cos^| < 0.8 to avoid 
any possibility of having results distorted by the extreme forward and backward regions 
which have low acceptance.^ The asymmetry is large, statistically significant and mass 
dependent. Within the framework of the partial wave analysis discussed below, the presence 
of only even values of £ would yield a symmetric distribution in cos 6*. Thus the observed 
asymmetry requires that odd-^ partial waves be present and that they interfere with even-£ 
partial waves to describe the data. Note that the decrease in asymmetry in the 1.4 GeV/c^ 
region can be (and will be shown to be) caused by the phase difference between the even-£ 
and odd-£ waves approaching 7r/2 rad. 



The azimuthal angular distribution as a function of tjti effective mass is shown in Fig. [15. 



The observed structure has a clear sin component which indicates the presence of m = 1 



'' The asymmetry function was plotted for various ranges of the decay angle and the presence of a strong 
asymmetry was noted in all cases. 



natural-parity-exchange waves in the production process. (See the discussion in Section [TVI 
below.) 

Shown in Fig.lTB| are the 7r~p and r]p effective-mass distributions for the data sample. 



It is important to note that the absence of baryon isobar production is required for the 
assumptions of our PWA to be valid. There is at most a very small amount of isobar 
production in the region M{7T~'p) < 2.0 GeV/c^ in Fig.p!6|a and none in Fig.|l6|b. The 
amplitude analysis described in Section |^ was checked to insure that isobar production did 
not effect our results. This was done by redoing that analysis after requiring M{tt^p) > 
2.0 GeV/c^ . The resulting intensities and phases did not change (other than an overall 
magnitude change due to the loss of events) in most cases by more than one standard 
deviation and in no case by more than 1.5 standard deviations. 

IV. PARTIAL- WAVE ANALYSIS 
A. Procedure 

A partial-wave analysis (PWA) |3lH33| based on the extended maximum likelihood 



method has been used to study the spin-parity structure of the t]tt^ system. We give in 
Appendices A and B some mathematical details regarding the techniques used in the partial- 
wave analysis. The formalism adopted in this analysis is somewhat different from those used 
by previous investigators. Although complete details used in the formalism are given in a 
recent publication by S. U. Chung [0, a portion of that work is reproduced in Appendices 



A and B in order to make this paper as complete and self-contained as possible. 

In Appendix A, a brief description of the formalism is given as are the relationships 
between the partial wave amplitudes (assuming i < 2) and the moments of the angular 
distribution. The technique of the extended maximum likelihood analysis is covered in 
Appendix B, where the interplay of the experimental moments, Eq. |Bl|, and the acceptance 
is described. (The experimental acceptance is incorporated into the PWA by using the 
accepted Monte Carlo events described above to calculate normalization integrals - see ref. 
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ii)- 

The partial waves are parameterized in terms of the quantum numbers J^*" as well as m, 
the absolute value of the angular momentum projection, and the reflectivity e [^. In our 



naming convention, a letter indicates the angular momentum of the partial wave in standard 
spectroscopic notation, while a subscript of means m = 0, e = —1, and a subscript of +(— ) 
means m = 1, e = +1(— 1). Thus, So denotes the partial wave having J^'-^Tif = 0"'"*'0^, while 
P_ signifies l""*"!", D^ means 2++1+, and so on. 

We consider partial waves with m < 1 in our analysis. This assumption is true in the 
limit of —t = 0, since the nucleon helicities give rise to the states with m = or m = ±1 
only. But this assumption can be dealt with — experimentally — since the moments H{LM) 
with M = 3 or M = 4 can be checked, to see how important the states \im) are in the data 
with \m\ > 2. This has been done with our data. The moments H{33), H{43) and if (44) 
are all small in the a^(1320) region, and a fit including |22) (shown in Section [iV (J| below) 
contains only a very small amount of this wave and is very broad. 

We also assume that the production spin-density matrix has rank one. This assumption 
is discussed in Appendix C. 

Goodness-of-fit is determined by calculation of a x^ from comparison of the experimental 
moments with those predicted by the results of the PWA fit. A systematic study has been 
performed to determine the effect on goodness-of-fit of adding and subtracting partial waves 
of J < 2 and M < 1. We find that although no significant structure is seen in the waves 
of negative reflectivity (see below), their presence in the PWA fit results in a significant 
improvement in goodness-of-fit compared to a fit which includes only the dominant positive- 
reflectivity partial waves. We have also performed fits including partial waves with J = 
3 and with J = 4. Contributions from these partial waves are found to be within one 
standard deviation of zero for most mass bins with M{riTT^) < 1.8 GeV/c^ and in all cases 
within two standard deviations of zero. Thus, PWA fits shown or referred to in this paper 
include all partial waves with J < 2 and m < 1 (i.e. So, Po, P-, Do, -D-, P+, and -D+). A 
non-interfering, isotropic background term of fixed magnitude determined as described by 
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Method 2 in Section lilDI is used. 



B. Results 

The results of the PWA fit of 38,200 events in the range 0.98 < M(r/7r-) < 1.82 QeV/c^ 
and 0.10 < |t| < 0.95 (GeV/c)^ are shown in Figs.|l^ and |l^. In Fig.|l^ the acceptance- 



corrected numbers of events predicted by the PWA fit for the i5+ and P+ intensities and the 
phase difference between these aniphtudes, A$, are shown as a function of M{riTx^). (The 
smooth curves shown in this figure are discussed below in Section [VB| .) There are eight 
ambiguous solutions in the fit []32| , p5| , |36| . These solutions are mathematically discrete but 
with equal likelihoods - that is, they correspond to exactly the same angular distributions. 
We show the range of fitted values for these ambiguous solutions in the vertical rectangular 
bar at each mass bin, and the maximum extent of their errors is shown as the error bar. 
These rectangular bars are quite small and thus not apparent for the Z)_|_ intensity, but they 
are quite clear for the P+-intensity and the phase-difference distributions. 



The a^(1320) is clearly observed in the D+ partial wave (Fig.p!7|a). A broad peak is seen 
in the P+ wave at about 1.4 GeV/c^ (Fig. |l^b). The phase difference A$ increases through 
the a^(1320) region, and then decreases above about 1.5 GeV/c^ (Fig.p!7|c). This phase 



behavior will allow us to study the nature of the P+ wave. (We note that there is a sign 
ambiguity in the phase difference and thus only the magnitude of A$ is actually measured.) 
Shown in Fig. |1^ are the fitted intensities for waves which are produced by negative- 
refiectivity (or unnatural-parity) exchange. The predicted numbers of events for these waves 
are generally small and are all consistent with zero above about 1.3 GeV/c^. Although there 
is some non-zero contribution from the D_ and (especially) the Sq waves below this region, 
the uncertainties and ambiguity ranges associated with these waves make it impossible to 
do a definitive study of them to determine their nature. In addition, the absense of a strong 
wave (such as the D+ wave in the natural-parity sector) to beat against these waves precludes 
us from drawing any conclusions about possible resonant behavior in the unnatural parity 
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sector. 

The forward-backward asymmetry noted earlier is due to interference in the natural- 
parity exchange sector rather than to the unnatural-parity exchange waves. This is illustrated 



in Fig.[T^ which shows the predicted asymmetry separately for the natural and unnatural 
parity exchange waves. It is clear that the asymmetry due to the unnatural-parity waves is 
about an order of magnitude less than that due to the natural-parity waves. Also shown in 



Fig. |T9| is the comparison of the asymmetry present in the data with that predicted by the 
fit. The fit clearly does an excellent job in representing the data points. 

The unnormalized spherical harmonic moments H{LM) and their prediction from the 
PWA fit as a function of mass are shown in Fig. ^ Here H{LM) = E^Ii Y^^Oi, (pi) {N being 
the number of events in a given bin of M{r]7r~)). The relationships between the moments 
and the amplitudes are given in Appendix A, Eqs. |A27| . All of the even moments (M = 0,2) 
are well-described by the fit as are most of the M = 1 moments. Some points with M = 1 
are somewhat less well fitted but, as we will discuss below, the significant conclusions which 
will be drawn from this work come from the natural-parity sector whose amplitudes do not 
contribute directly to the M = 1 moments (see Eqs. [A27| ). 

An examination of the if (30), if (32), ii(40) and ii(42) moments along with a comparison 
with Eqs. |A27| shows that the D+ amplitude dominates and demonstrates clearly that 
the P+ partial wave is required for the PWA fit to describe the experimental moments. 
These moments cannot be described solely by the combination of the D+ partial wave and 
experimental acceptance. 

The change in — log(Likelihood) as a function of the number of events for the P+ partial 
wave for the 1.30 < M{ri7i^) < 1.34 GeV/c^ bin is shown for all the ambiguous solutions in 



Fig. ^ The curves were obtained by fixing the P+ intensity at various values and maximizing 
the liklihood function varying all of the other parameters. All eight ambiguous solutions were 
found for each value of the P+ intensity. The number of predicted P+ events at the maximum 
of the liklihood ranges from 330 events to 530 events for the eight solutions with typical errors 
of 280 events. (A change in — log(Likelihood) of 0.5 corresponds to one standard deviation.) 
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For all solutions, the liklihood function gets so bad below 100 events that the P+ wave is 
clearly required to fit the data. Thus the observed variation in — log(Likelihood) further 
demonstrates that the P+ partial wave is required to describe our data. 

C. Systematic Studies 

PWA fits were performed for two different t ranges containing approximately equal num- 
bers of events. One bin spanned the range 0.10 < \t\ < 0.25 (GeV/c)^, and the other was for 
0.25 < \t\ < 0.95 (GeV/c)^. Both bins yielded comparable structures in the P_|_ wave, and 
the a^(1320) was the dominant feature of the D+ wave for both bins. The relative P+-D+ 



phase behavior for each bin was similar to the results for the integrated fit shown in Fig. pT| c. 

A PWA fit has been carried out excluding those events with \cos6gj\ > 0.8 (the region 
in which experimental acceptance is poorest). Neither the P+-wave intensity nor its phase 
variation relative to the D+ wave change by more than one standard deviation in any mass 
bin. 

As mentioned above, a PWA fit has been carried out including the natural parity exchange 
m = 2 amplitude (labelled D2+) in the fit. The results of this fit are shown in Fig.p2|. 
Comparing this fit with the fit shown in Fig. |l^ it is clear that the magnitude and phase 
behavior of the P+ wave is quite unaffected by inclusion of the m = 2 amplitude. 

To test for A and A^* contamination, a fit has been done in which events with M{tt~p) < 
2.0 GeV/c^ are excluded. As discussed earlier, the resulting intensities and phases did not 
change in most cases by more than one standard deviation and in no case by more than 1.5 
standard deviations. 

Fits were also carried out on Monte Carlo events generated with a pure D^ wave to 
determine whether structure in the P+ wave could be artificially induced by acceptance 
effects, resolution, or statistical fluctuations. Shown in Fig.^ are the results of such a fit. 
We do find that a P^ wave can be induced by such effects. This 'leakage' leads to a P+ wave 
that: (1) mimics the generated D+ intensity (and in our case would therefore have the shape 
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of the 02(1320)); and (2) has a phase difference A$ that is independent of mass. Neither 
property is present in our study and we conclude that the P+ structure which we observe is 
not due to 'leakage'. 

Fits have been performed allowing i = 3 and i = 4 waves. We find that these waves are 
negligible in the region below M{ri7T^) of about 1.7-1.8 GeV/c^, their intensities being less 
than one standard deviation from zero in almost all bins. (The largest number of events in 
any bin for the F_|_ wave was 34 ± 22 events and for the G+ wave was 140 ±100 events.) 

The data have been fit using different parametrizations of the background. The back- 
ground has been set at fixed values determined from the two different background estimates 
discussed previously. In another fit, the background has been set to zero. And finally, a fit 
was performed allowing the background level to be a free parameter. Although the nega- 
tive reflectivity waves do change somewhat for different treatments of the backgroundp^, the 
results for the -D+ and P+ waves and their relative phase do not change by more than one 
standard deviation in the entire region between 1.2 and 2.0 GeV/c^ except for a few isolated 
points which vary up to 1.5 standard deviations. 

D. Comparison with Previous Experiments 

These results for the P+ and D^ intensities and their phase difference are quite consistent 



with the VES results |jTj] as can be seen in Fig. ^ In particular, the behavior of the shape 
of the phase difference is virtually identical to that reported by VES.Q This is particularly 
noteworthy since, as will be seen below, it is this phase difference which allows us to draw 
conclusions regarding the nature of the P+ wave. 

Our results are compared with those of the KEK experiment |T^ in Fig.^. In this case. 



'^ Since the background and the ^o wave are both isotropic, the fitting program cannot distinguish 
between them. 

'^ The magnitude of the phase difference is shifted by about 20° relative to that of VES. A produc- 
tion phase shift would not be unexpected because of the differing energies and targets in the two 
experiments. 
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it is clear that the two resuhs differ. The KEK resuhs have a P-wave intensity which is 
narrower and a P-D phase difference which, within errors, is consistent with being constant 
as a function of M(?77r^). 

E. Comparison ^vith r] -^ tt+tt vr*^ Data Sample 

A second data set in another topological class (with two additional charged particles in 
the final state) has been used to study Reaction (|I|) with the decay mode rj -^ tt+tt^tt^. 
Besides having three forward charged particles instead of one, these events have a 7r° instead 
of an T] to be detected by the LGD. Since the vr" is one of three pions in the r] decay, its 
energy will be significantly less than the rj in the topology with only one charged track. Thus 
it is clear that the r] -^ tt+tt^tt^ data sample will have significantly different acceptance and 
systematics when compared to the rj ^ 2'-/ sample. 

Shown in Fig.|2^ is the tt+tt^tt'^ effective mass distribution from this data set. There is a 
clear r) peak as well as a strong peak in the u region. After kinematic fitting, a sample was 
obtained of 2,235 events which were consistent with Reaction |^ with rj —>■ tt+tt^tt'^. Fig.|27| 
shows the effective mass distribution for this sample of events. As expected, the a^(1320) 
dominates this spectrum. 

Although the data sample is very small for this topology, we have carried out an amplitude 
analysis in order to compare with the primary rj —* 2'-/ analysis. Results of the analysis are 



shown in Fig. ^ where we compare the shapes of the P+ intensities for the two data sets 
as well as the P+ — D^ phase differences. Despite the rather large statistical uncertainties, 
there is excellent agreement between these distributions. 

V. MASS-DEPENDENT FIT 

In an attempt to understand the nature of the P^ wave observed in our experiment, 
we have carried out a mass- dependent fit to the results of the mass- independent amplitude 
analysis. The fit has been carried out in the rin~ mass range from 1.1 to 1.6 GeV/c^. In this 
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fit, we have assumed that the D+-wave and the P+-wave decay amphtudes are resonant and 
have used relativistic Breit-Wigner forms for these amphtudes. 

A. Procedure 

We shall use a shorthand notation w to stand for the r/vr" mass, i.e. w = M{r]7i^). 
Representing the mass-dependent amplitudes for D+ and P+ as Ve{w) for i = 2 and 1, we 
may write 

Ve{w) = e'^' A,{w) B^) K + H^ - ^°) + c,{w - ^°)2]V' (2) 

where q is the r/vr" breakup momentum at mass w. Here (pi is the production phase (mass 
independent!^) associated with a wave i. The quantities Af(w) and -B^(q') are the standard 
relativistic Breit-Wigner form and the barrier factor, respectively, and are given below. The 
square-root factor has been introduced primarily to take into account possible deviations 
from the standard Breit-Wigner form, at values of w away from the resonance mass (denoted 
by wf). The overall normalization of a wave is governed by a^, while the constants hi and q 
allow for deviations in the mass spectra from the Breit-Wigner form. The constants a^, h^ 
and Q are all real, so that the square-root factor does not affect the rapidly varying phase 
implied by the standard Breit-Wigner form.Q 

The barrier functions ||3^, which are real, are given by: 



" We have tried a linear dependence in mass for the production phase; the fits did not require it. 

■'' For a Breit-Wigner form with a constant width, the phase rises 90 degrees over one full width centered 
at the resonance mass. 
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where z = {q/qj^)'^ and q^ = 0.1973 GeV/c corresponding to 1 fermi. Note that -B^(g) =^ 1 
as g ^ oo. 

The relativistic Breit-Wigner functions can be written 
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where F^ is the nominal width (mass independent) and Ti{w) is the mass-dependent width 
given by 
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where g° is the breakup momentum evaluated at w = w^. The mass- dependent phase shift 
SAw) is given by 
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and the overall phase for the £-wave amplitude is 



$^ = (t)i + 5Aw). 



We are dealing with two waves, P+ and -D+, and can only measure = 02 — 0i- Thus the 
phase difference being measured experimentally, corresponds to 

A^ = ^2 - ^1 = (f) + S2{w) - 6i{w) (9) 

Finally, the experimental mass distribution for each wave i is given by 

^ = \Veiw)\'pq (10) 

aw 

where pq is the phase-space factor for which p is the breakup momentum of the r/vr^ system 
(or of the final-state proton) in the overall center-of-mass frame in Reaction (|T]). Since the 
problem here is for a given ^/s, all other relevant factors, including that of the beam flux, 
have been absorbed into the amplitude itself, i.e. the constants a^, bg and q. 

The input quantities to the fit included, in each mass bin: the P+-wave intensity; the 
D+-wave intensity; and the phase difference A$ (the relevant formulas are given in @ and 
(|lO|)). Each of these quantities was taken with its error (including correlations) from the 
result of the amplitude analysis. One can view this fit as a test of the hypothesis that the 
correlation between the fitted P-wave intensity and its phase (as a function of mass) can be 
fit with a resonant Breit-Wigner amplitude. 

We find that the fit does not improve significantly when the P-f wave is modified from 
the Breit-Wigner form, and hence set 6i and Ci = for the final fit. We also note that the 
magnitudes of the quantities 62 and C2 in the final fit correspond to a small deviation of the 
D_|_-wave intensity of the order of 1%. 

B. Results 

Results of the fit are shown as the smooth curves in Fig. |l^a, b, and c. The mass and 
width of the J^^ = 2++ state (Fig. 0a) are (1317 ±1 ±2) MeV/c^ and (127 ±2 ±2) MeV/c^ 
respectively |3^. (The first error given is statistical and the second is systematic.) The mass 



and width of the J^^ = 1~+ state as shown in Fig. |T^ are (1370 ±16 +3J) MeV/c^ and (385 
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±40 J^^„g ) MeV/c^ respectively. Shown in Fig. p!7| d are the Breit-Wigner phase dependences 
for the 02^(1320) (hne 1) and the P+ waves (hne 2); the fitted D+— P+ production phase 
difference (hne 3); and the fitted D+— P+ phase difference (hne 4). (Line 4, which is identical 



to the fitted curve shown in Fig. p!7|c, is obtained as line 1 — line 2 + line 3.) 

The systematic errors have been determined from consideration of the range of solutions 
possible because of the ambiguous solutions in the PWA. Since there are 8 ambiguous solu- 
tions per mass bin and we are fitting over 12 mass bins, it is clearly impossible to try all 8^^ 
possible combinations. Instead, we have fit some 10^ combinations where the values to be 
fitted in each mass bin have been chosen at random from among the 8 ambiguous PWA solu- 
tions. The resulting fits generally clump into a group with reasonable values of x^ /dof{< 2) 
and into a group with poor values. The systematic errors on the mass and width given above 
are taken from the extremes observed for the solutions with reasonable values of 'x^ /dof . 
The central values quoted above are taken from a fit which uses the average values of the 
input parameters in each bin. 

The fit to the resonance hypothesis has a x^/dof of 1.49. The fact that the production 
phase difference can be fit by a mass-independent constant (of 0.6 rad) is consistent with 
Regge-pole phenomenologyp~| in the absence of final-state interactions. If one attempts to 
fit the data with a non-resonant (constant phase) P+ wave, and also postulates a Gaussian 
intensity distribution for the P+ wave, one obtains a very poor fit with a x^/dof of 7.08. 
Finally if one allows a mass- dependent production phase, a x^/dof of 1.55 is obtained for 
the non-resonant hypothesis — but in this case the production phase must have a very 
rapid variation with mass.[^ Furthermore, for this non-resonant hypothesis, one must also 
explain the correlated structure observed in the P+ intensity — a structure which is explained 
naturally by the resonance hypothesis. 



An attempt |3^ to explain our result as the interference of a non-resonant Deck-type 



^ The signature factor and the residue functions are at most t-dependent (not mass dependent) (see 
ref. 0). 

'* The fit requires a hnear production phase difference with a slope of -4.3 rad/GeV. 
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background and a resonance at 1.6 GeV/c^ can reproduce this correlation. (Evidence for an 
exotic meson with a mass near 1.6 GeV/c^ has been reported ^^ by our collaboration.) How- 



ever, this explanation is excluded because of the recent observation ||TB[ by the Crystal Barrel 
collaboration which confirms the presence of this state produced in nucleon-antinucleon an- 
nihilation. The Deck-effect is a mechanism applicable to peripheral production but not to 
annihilation. 

Our fitted parameters for the J^'" = 1""^ resonance are compared in Table |I| with the 
values reported by the Crystal Barrel experiment |T^. That experiment reports that a 
J^*^ = 1^+ resonance in the rjir channel is required to fit their data in the annihilation 
channel pn -^ Tr^vr^r/. Their fitted parameters are very consistent with those determined 
from our mass-dependent analysis. 

C. Other Systematic Studies 

1. Sensitivity to the D-wave Intensity Distribution Function 

In order to determine the sensitivity of the results of our mass- dependent analysis to 
the exact function being used to fit the D-wave intensity distribution we have redone the 
fit using two other hypotheses. First we have performed a fit in which the mass-dependent 
amplitude is given by Eq.(|^), but with 62 = C2 = 0. Second, we have taken 62 = C2 = in 
Eq.(^ and also replaced the Blatt-Weisskopf barrier functions for each wave by the factor 
q^. Although the resulting fits are poorer in quality, we find that the parameters of the fit 
do not change by amounts greater than the systematic uncertainty described above. 

2. Sensitivity to Leakage 

As shown in Section [iV C| , the P+ wave observed in our data is not consistent with 
"leakage" . That is, the analysis shows that the intensity and phase motion of the P+ wave 
do not have the characteristics of the wave which is artificially generated from a pure D^ 
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wave due to possible incomplete knowledge of the resolution or detection inefficiency. This 
does not preclude the possibility of some leakage being present in the data and distorting 
the results of the mass dependent analysis (MDA). In this section, we describe a test which 
has been carried out to study the sensitivity of our MDA results to possible residual leakage 
being present in the data. 

The fit which has been carried out is a mass dependent partial wave analysis (MDPWA). 
In such a fit, the PWA is carried out as in Section |^ but instead of carrying it out separately 
for each rjir^ mass bin, all bins are fit simultaneously and are tied together with a mass- 
dependent function for each partial wave. That is, the extended maximum likelihood function 
of the form given by Eq. [B16| is generalized to include mass dependence: 

n . 

ln£ oc Y. ln/(f^i, Wi) - / dQdw r/(fi, w) I{Q, w). (11) 

i 

The free parameters in the fit include, in addition to the amplitudes of the partial waves, 
the Breit-Wigner masses, widths, and intensities as well as mass-independent production 
amplitude phases. 

For simplicity, we have taken the Breit-Wigner form of Eq. ^ to describe each of the 
partial waves. A common mass and width were used for the D^, D_, and Dq partial waves. 
Similarly, the P+, P_, and Pq waves were assumed to be described by a common mass and 
width. The S wave was assumed to have its own mass and width. The constants bg and q 
of Eq. 1^ were taken to be zero for the Dq, D^, P_|_, Pq, P_, and 5*0 waves. Each wave was 
allowed to have its own normalization constant and production phase. 

In order to include leakage in the fit, a leakage amplitude Pi'^ with the characteristics 



obtained in the leakage study of Section |IVC| was defined. This amplitude was taken to have 
the shape of the D^ amplitude as well as its Breit Wigner phase dependence. Its production 
phase was fixed at (p^'"''^ = 80° This amplitude was then combined coherently with the P+ 
signal to give an effective amplitude given by the expression: 

Pi^"\w) = P4w)+Pf'\w) (12) 
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where 

p!^'\w) = af^ e*^'"' A,{w- w^, T^) B2{q)[l + bt{w- w^) +ht{w- wlfY'\ (13) 



The resuhs of the MDPWA fit are shown as the smooth curves in Fig. ^. Also shown 
as the points with error bars are the results of the mass-independent PWA. It is clear that 
the two analyses give consistent results. Shown in Fig. p9| a are the P'^ intensity (curve 
3) along with the P+ signal intensity (curve 1). The leakage is shown in curve 2 as the 
sum of the leakage intensity and the (signal - leakage) interference term. The fitted leakage 
contribution is equal to TZ^ = 0.018 where we have defined 

n, ^ [^) . (14) 

The results of the fit are given in Table ^ where they are compared with those of the 
combined PWA and mass dependent fit. The results are quite compatible when one takes 
into account the systematic errors. The biggest difference is in the fitted width of the P+ 
state which is larger for the MDPWA. 



In Fig. 30 are shown the fitted values for the mass (curve 1) and the width (curve 2) of 
the l""*" resonance as well as the change in \nC as a function of 7^^ (the leakage fraction). 
We note that the mass and width are very insensitive up to values of TZe = 5%, above which 
the fit becomes very unlikely. 

D. Cross Section Estimate 

In order to estimate the cross section for production of the observed P+ state (which we 
now refer to as 7ri(1400)), we have fitted published total cross sections pT|-[^ for 02^(1320) 



production to a function of the form a = A{p/pq) " + -B, where p is the beam momentum and 
Po is set to 1 GeV/c. Experiments with poor a]^(1260)/a2"(1320) separation were excluded 
from the fit. The best fit gave: A = 5099 ± 221 fib; n = 1.88 ± 0.03; and B = 39.2 ± 2.0 fib. 
From this we estimate the total cross section for 02^(1320) production at 18.2 GeV/c to be 
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61.1 ± 2.2 fib. This is in good agreement with the result of 62.56 ± 2.92 /x6 measured |^ at 
18.8 GeV/c. 

From the results of our PWA, we find that, in the rj mass range from 1.10 to 1.58 GeV/c^ 
there are 60, 332 ± 2, 060 D+ events, and there are 3, 321 ± 1, 245 P+ events. Here the error 
for the number of D^ events is statistical and the error for the number of P+ events includes 
uncertainties due to ambiguities. One thus obtains: cr{n^p -^ pvrf (1400)) *i?i?(7rf(1400) -^ 
riTc-) = 0.49 ± 0.19 /x6. 

VI. SUMMARY AND CONCLUSIONS 

In this paper, we have discussed the details of the amplitude analysis of data from Reac- 
tion 0. Interference between D-wave and P-wave amplitudes produced with natural parity 
exchange is required in order to explain the data. Using this interference, we have shown 
that the P-wave phase has a rapid variation with mass and that this phase variation coupled 
with the fitted P-wave intensity distribution is well-fitted by a Breit-Wigner resonance with 
mass and width of (1370 ±16 +^J) MeV/c^ and (385 ±40 +^^^5) MeV/c^ respectively. Since 
a P-wave resonance in the rjTT system has J^*" = 1 ^, it is manifestly exotic. The exact 
nature of the observed state awaits further experimentation. 
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APPENDIX A: PARTIAL- WAVE FORMULAS 

In this appendix, the angular distributions are derived for the rjir^ system produced in 
Reaction ([1|). The distributions are given both in terms of the moments and the amphtudes 
in the reflectivity basis. For a system consisting of S, P and D waves, explicit formulas for 
the moments as functions of the partial waves are also given. 

In the Gottfried- Jackson (GJ) frame, the amplitudes may be expanded in terms of the 
partial waves for the r^vr" system: 

f/fc(fi) = ^ VemkAimin) (Al) 

where Vgmk stands for the production amplitude for a state \im) and k represents the spin 
degrees of freedom for the initial and final nucleons {k = 1,2 for spin-nonfiip and spin- flip 
amplitudes). Aim{^) is the decay amplitude given by 



AU^) = ^^ Di*(0, 9, 0) = Yrm (A2) 

where the angles Q = {6, 0) describe the direction of the rj in the GJ frame. It is noted, 
in passing, that the small (i-function implicit in (|A2|) is related to the associated Legendre 
polynomial via 



dUo) = {-y 



The angular distribution is given by 



^^ "^^^ Pr(cos9) (A3) 



\(^ + 



m] 



iin) = Y,\Ukm' (A4) 

k 

It should be emphasized that the nucleon helicities are external entities and the summation 
on k is applied to the absolute square of the amplitudes. A complete study of the rjn' 
system requires four variables: M(?77r~), —t and the two angles in Q. The distribution 
is therefore to be applied to a given bin of M{riTT^) and of — t. 
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The angular distribution may be expanded in terms of the moments H{LM) via 

n^) = E C-^) H{LM)D^^*{<l>, e, 0) (A5) 



where 



H{LM) = Y: -^^ pZm'{i'm'LM\im){i'OLOm (A6) 

im \ I 



i'm' 



where p is the spin-density matrix given by 

Pmm' = E ^irnkVi*m'k (A7) 

k 

It is seen that the moments H{LM) are measurable quantities since 

H{LM) = Jdni{n)Dl,^{(l),e,o) (as) 

The normalization integral is 

H{00) = f dVt I{Q) (A9) 

The symmetry relations for the /J's are well-known. From the hermiticity of p, one gets 

H*{LM) = (-)*^ H{L - M) (AlO) 

and, from parity conservation in the production process, one finds 

H{LM) = {-f^ H{L - M) (All) 

These show that the if 's are real. 

Parity conservation in the production process can be treated with the reflection operator 
which preserves all the relevant momenta in the S'-matrix and act directly on the rest states 
of the particles involved. It is important to remember that the coordinate system is always 
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defined with the y-axis along the production normal. In this case the reflection operator is 
simply the parity operator followed by a rotation by vr around the y-axis. 
The eigenstates of this reflection operator are 

\eim) = e(m){\im) - e(-)"^|^ - m)} (A12) 



where 
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(A13) 



For positive reflectivity, the 171 = states are not allowed, i.e. 

leiO) =0, if e = + (A14) 

The reflectivity quantum number e has been defined so that it coincides with the naturality 
of the exchanged particle in Reaction ([l|). One can prove this by noting that the meson 
production vertex is in reality a time-reversed process in which a state of arbitrary spin- 
parity decays into a pion (the beam) and a particle of a given naturality (the exchanged 
particle) 

JVj ^ s'7. + TT (A15) 



where r/'s stand for intrinsic parities. The helicity-coupling amplitude F'^ for this decay |^ 
is 

A'^{M)^F^Di;,{<p,,e,,0) (A16) 

where A is the helicity of the exchanged particle and the subscript p stands for the 'produc- 
tion' variables. M is the z-component of spin J in a J rest frame. From parity conservation 
in the decay, one finds 
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Fi = -Fl, (A17) 

where one has used the relationships rjj = (— )'^ (true for two-pseudoscalar systems) and 
T]s = (— )* (natural-parity exchange). The formula shows that the helicity-coupling amplitude 
F"^ is zero if A is zero. Since angular momentum is conserved, its decay into two spinless 
particles in the final state cannot have M = along the beam direction (the GJ rest system), 
i.e. the D"^-function is zero unless M = X, ii 9p = (j)p = 0. Finally, one may identify J with 
i and M with m, which proves (|A14|) . 



The modified Z)-f unctions in the reflectivity basis are given by 

'Di:,{<f>, e, 0) = e{m) [DZ{<f>, e, o) - e{-r DiU<j^, e, o)] (ais) 

It is seen that the modified D-functions are real if e = — 1 and imaginary if e = +1: 

^-^DZi<P,0,O) = 2e{m)dUO)cosm<p 
Wd^*(0, e, 0) = 2t9im)dU9) sinm0 



(A19) 



The overall amplitude in the reflectivity basis is now 

'Uk{n) = J2 'Vimk 'A,^{Q) (A20) 



Irn 



where 



'AeU^) = f-^'Dlt{<P,9,0) (A21) 

and the resulting angular distribution is 

/(fi)=5:rf/fc(fi)|2 (A22) 

It is seen that the sum involves four non-interfering terms for e = ± and k = 1,2. The 
absence of the interfering terms of different reflectivities is a direct consequence of parity 
conservation in the production process. We use the partial wave amplitude notation 
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[^]o = ^"V,o, [£]- = (-Vn, [£]+ = W\/a (A23) 

where [i] stands for the partial waves S, P, D, F and G for i =0, 1, 2, 3 and 4. 

Consider an example where the maximum £ is 2. One sees that there are in general twelve 
possible non-zero experimental moments: 



H{00), H{10), H{11), H{20), H{21), H{22) 
H{30), H{31), H{32), H{40), H{41), H{42) 



(A24) 



while the partial waves [i] are, for unnatural-parity exchange, 

So, Po, P-, Do, D_ (A25) 

and, for natural-parity exchange, 

P+, D^ (A26) 

One wave in each naturality can be set be real (5*0 and P+, for example), so that there are 
again twelve real parameters (to be determined). It is helpful to write down the moments 
explicitly in terms of the partial waves: 
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H{10) = -^SoPo + -^PoDo + -^{P-D- + P+D+) 

Hi20) = ^SoDo + Ip^ - liP'_ + Pi) + ^-Dl + \{Dl + Dl) 

H{22) = \^-{Pl-Pl) + \J\{Dl-Dl) 

3 (A27) 

^(30) = — =(v^Po^o - P-D. - P+D+) 

H{31) = iy'^(2Po^- + ViP^Do) 

HiAl) = \\JIdoD. 
H{A2) = y^{Dl - Dl 



^l) 



APPENDIX B: MAXIMUM-LIKELIHOOD ANALYSIS 

This appendix is devoted to an exposition of the experimental moments, the acceptance 
moments and the acceptance-corrected (or 'true') moments and the relationships among 
them. Finally, the extended likelihood functions are given as functions of the 'true' and 
acceptance moments. 

One may determine directly the experimental moments (unnormalized) as follows: 

n 

H,{LM) = Y^ Di,,{<P,,,e,,0) (Bl) 

i 

where the sum is over a given number n of experimental events in a mass bin. But this is 
given by, from ( |A8| ) , 
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H^{LM) = Jdn r]{n) i{n) Di,^{<p, e, o) (B2) 

where r]{fl) represents the finite acceptance of the apparatus, and includes software cuts, if 
any. From ([A5[ ), one finds that 

H^{LM) = J2 H{L'M') ^,{LML'M') (B3) 

L'M' 

where 

/or' + i\ r 
^^{LML'M') = ( ^^1 J dfi Vm D^,,,i<p,e,0)D^,',*{<p,e,0) (B4) 

Note that the \l/'s have a simple normalization 

^.,{LML'M') = 5ll'Smm' (B5) 



in the limit ri(il) = 1. The integral ( [B4| ) can be calculated using a sample of 'accepted' MC 
events. Let A^^^ be the number of accepted MC events, out of a total of A^ generated MC 
events. Then, the integral is 

/ 2 /' ' + 1 \ 1 ^"^ 
^^{LML'M') = (^^1 ^E DiU<p,,9,,0)D^^*{cP,,9,,0) (B6) 

Equation ( |B3| ) shows that one can predict the experimentally measurable moments (|Bl|), 
given a set of true moments {H} and the \E''s; this provides one a means of assessing the 
goodness of fit by forming a x^ based on the set {H^}. 

There exists an alternative method of determining \E''s. For the purpose, one expands 
the acceptance function t]{Q) in terms of the ortho normal D-functions, as follows: 

V{n) = E(2^ + ^nLM) Di,t (0, e, 0) (B7) 

LM 

where ^(LM) is given by 

^LM) = ^Jdn 7^{n) D^,o(0, 0, 0) (B8) 
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The complex conjugate is, from the defining formula above, 

C{LM) = {-Y^i{L -M) (B9) 

so that the acceptance function can be made explicitly real 

r^iSl) = Y,(2L + l)r(M)Re [i{LM) Di,*{4>, 9, 0)} (BIO) 

LAI 

where 

t{M) = 2, M > 0, 

= 1, M = 0, (Bll) 

= 0, M < 



One sees that t{M) = 46*2 (M) where e{M) is defined in Eq. i \KT^ 



A set of ^{LM) specifies completely the acceptance in the problem. The normalization 
for the acceptance function has been chosen so that a perfect acceptance is given by 77 (fi) = 1 
and ^(LM) = 6loSmo- The ^(LM)'s can be measured experimentally using the accepted MC 
events 

aLM) = ^^ E ^Mo(0., 0^, 0) (B12) 



Finally, substituting ( p7|) into ( p^ , one finds 

<^^{LML'M')= Y. (21" + 1)C{L"M"){LML"M"\L'M'){L0L"0\L'0) (B13) 

L"M" 

This formula shows an important aspect of the ^{LM) technique of representing acceptance. 
Although (p^) involves a sum in which L and M could be extended to infinity for an arbitrary 
acceptance, there is a cutoff if the set {H} has maxima Lm and M^ [see (p3|) ]. The formula 
above demonstrates that L" < 2Lm and \M"\ < 2Mm- 

In a partial-wave analysis, it is usually best to take a set of the partial waves, [£]o, [C\~ 
and [£] + , as unknown parameters to be determined in an extended maximum-likelihood fit. 
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Since there is an absolute scale in an extended maximum-likelihood fit, one then has the 
predicted numbers of events for all the partial waves, corrected for finite acceptance and 
angular distributions. The partial waves in turn give rise to a set of predicted moments 
{H}. But the H{00) is not 1 but the total predicted number of events from the fit, i.e. one 
should be using the unnormalized moments. It is possible to choose iJ's as unknowns in the 
fit, but the two sets of iJ's should be the same ideally — this affords one an effective way of 
assessing self-consistency between the chosen moments and the partial waves. 

For completeness, a short comment is given about the extended likelihood functions. 
The likelihood function for finding 'n' events in a given bin with a finite acceptance 1]{Q) is 
defined as a product of the probabilities. 



C oc 



n" 



n\ 



n 



ji{n)ri{n) dn 



(B14) 



where the first bracket is the Poisson probability for 'n' events. This is the so-called extended 
likelihood function, in the sense that the Poisson distribution for 'n' itself is included in the 
likelihood function. Note that the expectation value n for n is given by 



noc f I{n)r]{n) dVL 



(B15) 



The likelihood function C can now be written, dropping the factors depending on n alone. 



C oc 



n Ii^^ 



exp 



i{n) 7]{vl) dn 



The 'log' of the likelihood function now has the form, 

n „ 

ln£ oc J2 ln/(fii) - dQ r]{Q) I{Q) 



(B16) 



which can be recast in terms of the ^(LM)'s 

n 

ln£ oc ^ ln/(a) - ^ (2L + 1) H{LM) C{LM) 

i LM 

n 

OC ^ ln/(a) - ^ (2L + 1) t{M) H{LM) Re^(LM) 



(BIT) 



LM 
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H{LM)^s may be used directly as parameters in the fit or may be given as functions of tfie 
partial waves. It is interesting to note that the ^(LM)'s for L > L^ and |M| > M^ are not 
needed in the likelihood fit. Note also that only the real parts of the ^(LM)'s are used in 
the fit. 

It should be borne in mind that a set of the moments {H} may not always be expressed 
in terms of the partial waves. This is clear if one examines the formulas ( |A271 ). Consider, 
for example, an angular distribution in which H{10) is the only non-zero moment. But this 
moment is given by a set of interference terms involving even-odd partial waves. So at least 
one term cannot be zero — for example, the interference term involving S- and P-waves. But 
then neither if (00) nor if (20) can be zero, since both S- and P-waves are non-zero. One 
must conclude then that a x^ based on the set {H^} may not necessarily be zero identically. 

APPENDIX C: RANK OF THE DENSITY MATRIX 

An assumption needed for the partial- wave analysis is that the density matrix has rank 1, 
i.e. the spin amplitudes do not depend on the nucleon helicities. Our justification, so far, has 
been that the fitted partial waves are very reasonable, that these waves can be fitted with 
a very simple mass- dependent formula, that Pomeron-exchange amplitudes are in general 
independent of nucleon helicities, and so on. . . . 

The purpose of this appendix is to point out that, under a simple model for mass depen- 
dence of the partial waves, it is possible to prove that the spin density matrix has rank 1. 
Suppose that one has found a satisfactory fit under a rank-1 assumption. One can then show 
that, even if the problem involves both spin-nonflip and spin-flip at the nucleon vertex — i.e. 
it appears to be a rank-2 problem — the spin density matrix in reality has a rank of 1. Al- 
though this note is based on the results of our r^vr^ analysis, the derivation does not depend 
on the decay channels; the conclusions apply equally well to other decay channels. 

This note relies on some technicalities generally well known, and so they have been 
presented without attribution. The reader may wish to consult a number of preprints and/or 
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papers, which deal with them in some detail pT|j3^,|52|-p4 



1. Partial Waves Produced via Natural-parity Exchange 

Consider the rj7r~ system produced via natural-parity exchange. It consists mainly of just 
two waves -D+ and P+ in the a^(1320) region. Assume these are the only waves. Without 
loss of generality, the decay amplitudes [^ can be considered real, i.e. 



A^{n) = \ —V^dl^i^d) sin0= -J— V3 sin0cos^sin0 
V 47r V 47r 

A^{n) = J^ V2dlf^{9) sin0 = -J^ sin^ sin0 



(CI) 



Since one deals with the partial waves produced only by natural-parity exchange, one can 
drop the subscript '+' from the waves, and the angular distribution is simply given by 



3 \ ^_ „ _ 2 

Air 
3 



oc 



oc 



An 



yhD cos 6 + P sin^ 9 sin^ 

'SIZ^I^ cos'^e + 2V5^{D*P} cos^+|P|2 



(C2) 



sin 9 sin^ 



The integration over the angles can be carried out easily, to obtain 



f I{n)dn(x \D\^ + \P\ 



(C3) 



as expected. 

It is easy to calculate the forward-backward asymmetry A{F, B) = {F — B)/{F + B 
(see section |IlD 

3^5 ,._ \P\\D 



A{F,B) = ^— cos(A$) 



4 



(|P|2 + |D|2) 



(C4) 



where A<l> is the phase difference between the P and D waves. 
The spin density matrix is given by 

I{Q) (X \DA^{Q) + PA^{Q)\^ = ^ Pk,k'AkAl, 

k,k' 



(C5) 
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where {k, k'} = {1, 2} and '1' ('2') corresponds to D (P). From this definition, one sees that 

/ |D|2 DP*\ 

P=[ , J (C6) 

\D*P |P|2 / 

One can work out the eigenvalues of this 2x2 matrix: 

A = {|Dp + |Pp, 0} (C7) 

One of the two allowed eigenvalues is zero, i.e. the rank of this matrix is 1. This is the 
'rank-1' assumption one makes to carry out the partial- wave analysis and is valid for a given 
mass bin. 

Suppose now that the rank is 2, i.e. 

i{n)oc\D,A^{n) + p,A^{n)\^ + \D2A^{n) + P2A^{n)\^ (cs) 

where subscripts 1 and 2 stand for spin-nonfiip and spin-flip amplitudes at the nucleon vertex 
for Reaction ([l|). Comparing (|C2|) and ( |C8|) , one flnds immediately 



|P|2=|p,|2 + |p2|2 (C9) 

^{P* D} = 3fJ{Pi* D^} + ^{P; D^} 

Let w be the effective mass of the r]7T~ system. If the mass dependence is included 
explicitly in the formula, one should write, in the case of rank 1, 

^^^ oc \D{w) A,{Q) + P{w) A,{Q)f pq (CIO) 

where p is the breakup momentum of the "qn^ system in the overall CM system and q is the 
breakup momentum of the rj in the rjir' rest frame. Note that both p and q depend on w. 
Note also that the w dependences of the partial waves D and P are given in the formula. 
Obviously, a similar expression could be written down for the case of rank 2. 
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One is now ready to make the one crucial assumption for a mass-dependent analysis of 
the D and P waves: one assumes that two resonances — in D and P waves, respectively — are 
produced in both spin-nonflip and spin-flip amplitudes. One may then write, for the rank-1 
case, 

D{w) = ae'"e'^^ sin^a 

(Cll) 

P{w) = he'^'' sm5b 

where a, h and the production phase a are all real and independent of the 777?" mass. In 
addition, one can set a > and 6 > without loss of generality. Here 5a and 5}, are the 
phase-shifts corresponding to the resonances and are highly mass dependent. In its generic 
form, the Breit-Wigner formula is given by the usual expression 

cot5 = ^i^ (C12) 

where Wq and Fq are the standard resonance parameters. In this note, the width is considered 
independent of w. Likewise, the barrier factor dependence for D and P is ignored.^ 
The formulas (|C11| ) are generalized to the case of rank 2, as follows: 

Diiw) =aie^"^e^'''^ Sin5a 

P^{w) =6ie*^' sin^fe 

(C13) 

D2H =a2e*"^e*'''" sin^^ 
P^{w)=b2e'^'- sin 6b 

Once again, a^, 6j and a^ are real, a^ > and hi > 0, and independent oi w. One finds, using 



Although simphfied formulas are used in this note, the results given here do not change even when 
correct formulas are used. Note that, to go over to a correct formulation for each wave, one needs to 
substitute the absolute value of the Breit-Wigner formula as follows: 



sin(5(it;) -^ B{q) 



T{w) 



sm5{'w) 



where B{q) is the barrier factor and V{w) is the mass-dependent width. It should be noted that the 
correction factors are all real, by definition. 
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2 2 2 

a = ai + a2 

b-' = bl + hi (C14) 

a6cos(a + 5a- Sb) = cti&i cos(q;i + 6a - Sb) + 0262 cos(a2 + Sa - Sb) 

A plot of cos(a + Sa — 6b) as a function of w is shown in Fig. ^ for three values of a, i.e. 
0°, 45° and 90°. The resonance parameters for a and b of 1.0 and 0.151 are taken from the 
mass dependent fit of Section |V B| as are the resonant masses and widths.^ The normalized 



absolute squares of the Breit-Wigner forms are given in Fig. 32, as is the 'normalized' inter- 



ference term. The same quantities, as they appear in Ref. [|I[, are shown in Fig. ^. This 
figure shows how important the interference term is compared to the P-wave term. Note 
also how rapidly the interference term varies as a function of w in the 02 (1320) region. This 
term, of course, is intimately related to the asymmetry in the Jackson angle and vanishes 



when integrated over the angle, i.e. it does not contribute to the mass spectrum [see (|C^ ) 



and (P3|)]. Fig. |3^ shows the contour plot of the intensity distribution in w vs. cos 6'; note 
the variation of the asymmetry as a function w. 

For the last equation in ( |C14| ) to be true for any mass, the coefficient of cos{6a — Sb) or 



sin(5a — 6b) on the left-hand side must be equal to that on the right-hand side, so that 

ab cos a = ai&i cos ai + 02^2 cos 02 

(C15) 

ab sin a = aibi sin ai + 0262 sin 0^2 

Taking the sum of the squares of the two formulas above and introducing the first two 
equations of (|C14| ), one obtains: 



2ai6ia2&2 cos ai cos 0:2 + 2ai6ia2&2 sin ai sin 02 

--ajbl + ajbj (C16) 

=0^62 (cos^ ai + sin^ ai) + a\b\ (cos^ 0^2 + sin^ 0^2) 



^ The value of a as determined from this fit is 37.46°; for the purpose of illustration, one may consider 
a = 45° close enough. 
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which is recast into 

= (0162 cosai — 0261 cos 0:2)^ + (ai&2 sinai — 0261 sin 02)^ (C'17) 



It is clear that each term must be set to zero, so that 

ai\ fO'2\ 

-— cosai =1— cos 02 

biJ V 62 / 

ai\ . fa2\ . 

-— sm «! = -— sm ^2 
hi J V 02 / 



(C18) 



Placing these back into (|C15| ), one deduces that 



a\ fai\ (a2 

— I cos ct = I 7- ) cos c^i = I 7- ) cos 0:2 

a\ . /ctA . (0.2 



(C19) 



One may take — alternately — the sum of the squares of the two formulas above, or a division 
of the second by the first, and obtain (remembering that the a's and 6's are non-negative 

real quantities), 

a Oi 02 

h h &2 (C20) 

tan a = tan ai = tan 02 

The last equation above demands that ai and ^2 are determined (up to ±7r), but they have 
to satisfy (|C19|) . It is therefore clear that one must set a = ai = 02. Next, one introduces 
two new real variables a; > and y > 0, given by 

(C21) 



X 

a b 

02 62 

a 



with the constraint x"^ + y"^ = 1. 

Now one can prove that the case of rank 2 is reduced to that of rank 1. Indeed, one sees 
immediately that 



and ( |C8| ) becomes identical to ( P2| 
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2. Discussion 

It is shown in this appendix that the problem of two resonances in D^ and P+ in the r/vr" 
system in (|1|) is — effectively — a rank-1 problem. For this to be true, the following conditions 
have to be met: 

(a) There exist two distinct resonances with different masses and/or widths. Note that the 
crucial step, from (|C14|) to (|C15| ), depends on that fact that 5a — 5f, is non-zero and is 



mass dependent. 

(b) There exists a satisfactory rank-1 fit with two resonances in a given mass region, in 
which each amplitude for D+ or P+ has the following general form 

7Wfc(w;,fi) = rfce*°^e^^'=("')/,HAfe(f]) (C23) 

where k = {1,2} and '1' ('2') corresponds to D^ (-P+)- 5k{w) is the Breit-Wigner 
phase and highly mass dependent, while r^ and a^ are mass independent in the fit. 
Of course, one of the two a^'s can be set to zero without loss of generality, so that 
there are three independent parameters, e.g. ri, r2 and ai (these were denoted a, b 
and a, respectively, in the previous section). fk{w) contains the absolute value of the 
Breit-Wigner form, plus any other mass- dependent factors introduced in the model. 
y4fc(n) carries the information about the rotational property of a partial wave k. 

(c) The same two -D+ and P+ resonances are produced in both spin-nonflip and spin-flip 
amplitudes, with the same general form as given above — but with arbitrary r^'s and 
afc's for each spin-nonflip and spin-flip amplitude. It has been shown in this appendix 
that only one set of r^'s and a^'s, i.e. ri, r2 and ai, is required for both spin-nonflip and 
spin-flip amplitudes. (This is indeed a remarkable result; the rank-2 problem entails a 
set of six parameters, but it has been shown that the set is reduced to that consisting 
of just three.) Therefore, the distribution function in both w and Q is given by 



40 



^-^^l^M.(u.,a)\'p, (C24) 



independent of the nucleon helicities. 

In another words, the spin density matrix has rank 1. The key ingredients for this result 
are that both spin-nonflip and spin-flip amplitudes harbor two resonances in D^ and P+ and 
that the production phase is mass-independent. It should be emphasized that the derivation 
given in this note does not depend on the existence of a good mass fit; it merely states that 
any fit with a mass-independent production phase is necessarily a rank-1 fit. Of course, the 
point is moot, if there exists no satisfactory fit in this model. 
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FIGURES 
FIG. 1. Experimental layout for E852. The nomenclature is defined in the text. 

FIG. 2. Average acceptance vs. riTT~ effective mass. 

FIG. 3. Average acceptance vs. cos6 for different ryvr" effective mass regions. 

FIG. 4. Average acceptance vs. (p for different r/vr" effective mass regions. 

FIG. 5. Average acceptance vs. \t'\ integrated over all r/vr" effective masses. 

FIG. 6. Two-photon effective mass distribution for events in the a^(1320) effective-mass region. 
The central cross-hatched region shows the events which remain after SQUAW fitting. The shaded 
sidebands show the regions selected to estimate the background using Method 1 (see text). 

FIG. 7. Missing-massed squared distribution. The dashed histogram shows the distribution of 
events which remain after kinematic fitting which leads to a rather sharp cutoff. 

FIG. 8. Missing-mass squared vs. two-photon mass. The elliptical regions are used to estimate 
the background using Method 2 (see text). 

FIG. 9. Effective-mass distribution of the background estimated from the rj sidebands (Method 
!)• 

FIG. 10. Angular distribution of the background (Method 1) shown separately for (a.) the 
low-mass sideband and (b.) the high-mass sideband of the rj. Events are plotted which fall in the 
a2^(1320) effective-mass region. 

FIG. 11. rjTT^ effective mass distribution uncorrected for acceptance. The shaded region is an 
estimate of the background using Method 2. 
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FIG. 12. Distribution of \t'\ = \t — tminl (acceptance-corrected). This experiment (solid dots) 
compared to a second experiment (open circles, see text) and to a Regge Pole fit. 

FIG. 13. Distributions of the acceptance-corrected cosine of the decay angle in the GJ frame 
for various effective mass selections. 

FIG. 14. Forward-backward asymmetry (acceptance-corrected) as a function of effective mass. 
The asymmetry = {F — B)/{F + B) where F(B) is the number of events for which the rj decays in 
the forward (backward) hemisphere in the GJ frame. 

FIG. 15. Distributions of the acceptance-corrected Treiman-Yang angle (j) in the GJ frame for 
various effective mass selections. 

FIG. 16. Effective mass distributions for: a.) the vr^p; and b.) the r]p systems for the final 
event sample (uncorrected). 

FIG. 17. Results of the partial wave amplitude analysis. Shown are a) the fitted intensity 
distributions for the D^ and b) the P+ partial waves, and c) their phase difference A<I> . The 
range of values for the eight ambiguous solutions is shown by the central bar and the extent of the 
maximum error is shown by the error bars. Also shown as curves in a), b), and c) are the results of 
the mass dependent analysis described in the text. The lines in d) correspond to (1) the fitted D^ 
Breit-Wigner phase, (2) the fitted P+ Breit-Wigner phase, (3) the fitted relative production phase 
(j), and (4) the overall phase difference A$. 

FIG. 18. Results of the partial wave amplitude analysis. Shown are the fitted intensity distri- 
butions for the waves produced by unnatural-parity exchange. 

FIG. 19. Forward-backward asymmetry as a function of effective mass. Shown are: the total 
asymmetry in the data (closed circles); the predicted asymmetry from the PWA fit (open squares); 
the prediction of the fit for that part of the asymmetry due to natural-parity exchange (filled 
squares); and the prediction of the fit for that part of the asymmetry due to the unnatural-parity 
exchange waves (open circles). 
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FIG. 20. Experimental moments H{L, M) (open circles) shown with the predicted moments 
(open triangles) from the amplitude analysis. 

FIG. 21. Value of the log likelihood as a function of the number of P^ events in the PWA fit for 
all 8 ambiguous solutions. The inset shows a view with expanded scales. Because some solutions 
are very close to each other, not all 8 solutions are distinguishable on this figure. 

FIG. 22. Results of the partial wave amplitude analysis when the natural parity exchange m = 2 
amplitude is included. Shown are the fitted intensity distributions for a) the D+, b) the ?+ , and 
d) the D2+ partial waves. Shown in c) is the phase difference A(/> between the 0+ and the P+ 
partial waves. 

FIG. 23. Fitted P+ intensity and P+—D^ phase difference for the Monte Carlo sample generated 
with a pure D+ sample of a^(1320) events. 

FIG. 24. Comparison of the results of this amplitude analysis with the VES experiment. Shown 
are the P+ — I?+ phase difference and the P+ intensity as a function of r]TT~ effective mass for each 
experiment. Note that the left-hand scales are for E852 and the right-hand scales are for VES. 

FIG. 25. Comparison of the results of this amplitude analysis with those of the KEK experiment. 
Shown are the P+ — Dj^ phase difference and the P+ intensity as a function of r/vr" effective mass 
for each experiment. Note that the left-hand scales are for E852 and the right-hand scales are for 
KEK. 

FIG. 26. The vr+Tr^vr'' effective mass distribution for events with the topology of three forward 
charged tracks, one recoil charged track, and two photon clusters consistent with a vr*^. 

FIG. 27. The r/vr^ effective mass distribution for the 77 -^ vr+vr^Tr*^ event sample. 

FIG. 28. Comparison of the results of the amplitude analysis for the 77 -^ vr+Tr^vr*^ (filled 
triangles) and the r/ ^ 27 (open circles) samples. The ordinate scale for the P+ intensity is for the 
?7 ^ 27 fit only. Thus, only the shapes of the P+ intensity distributions should be compared. 
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FIG. 29. The fit results of the MDPWA (curves) and of one solution for the mass indepen- 
dent PWA (crosses) for the r]TT^ system: a) P+ , b) D-^- intensities and c) their relative phase 
A(/)(P_|- — D+). Fig. la also shows the contributions of the 1 ^ signal intensity (1), the sum of the 
leakage and (signal - leakage) interference term (2) and the complete 1 ^ wave (3). 

FIG. 30. Dependence of the 1 ^ signal parameters Wi (1), F^ (2) and the change in the 
-log(Likelihood) function relative to its minimum (3) on the leakage contribution TZi at </> = 80°. 
The black points are at the position of the likelihood extremum. 

FIG. 31. cos(a + 6a - 4) as a function of w from 1.2 to 1.6 GeV/c^ for a = 0°(o), a = 45°(+) 
and a = 90°(n). 

FIG. 32. sin^ 5a (o), sin^ Sb (+) and sin (5a sinJf, cos(a + Sa — S^) (□) as a function of w from 1.2 
to 1.6 GeV/c^, using a = 45°. 

FIG. 33. a^ sin 6a (o), b'^ sin 6b (+) and 2ab sin (5a sin6b cos(q + 6a — 6b) (n) as a function of 
w from 1.2 to 1.6 GeV/c^, where one has assumed that a = 1.0, b = 0.20 and a = 45°. 

FIG. 34. Angular distribution in cos^ as a function of w from 1.2 to 1.6 GeV/c^, where one 
has assumed that a = 1.0, b = 0.151 and a = 37.46°. 
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TABLES 



TABLE L Reduction in the data sample as a function of the cut type. 



Cut 


Remaining number of events 


Fraction removed (%) 


Number of triggers 


47 X 10^ 


- 


Topological and trigger cuts 


583,094 


98.8 


r] preselection (C.L. > 10""^) 


270,364 


53.6 


Removal of runs with LGD 


159,871 


40.9 


trigger processor failure 






LGD fiducial cut 


146,584 


8.3 


Photon-hadron distance cut 


145,710 


0.60 


Missing-mass-squared cut 


103,341 


29.1 


Confidence level cut 


85,888 


16.9 


A(/.<8° 


53,219 


38.0 


TPX2 cut 


49,113 


7.7 


Cut on position at EV/BV 


47,235 


3.8 


0.10 < \t\ < 0.95 (GeV/c)2 


42,676 


9.7 


0.98 < M{r]7r-) < 1.82 GeV/c^ 


38,272 


10.3 



TABLE IL Comparison of the results of E852 and the Crystal Barrel for the parameters of the 
J^'^ = 1 ^ resonance. 



Mass (MeV/c^ 



Width (MeV/c2 



E852 

Crystal Barrel 



1370 ±16 tlo 



1400 ±20 ±20 



385 ±40 



310 ±50 



+65 
-105 

+50 
-30 
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TABLE III. Comparison of the results of the PWA combined with a separate mass dependent 
fit (MDF) with those of the MDPWA with leakage. 



Meson 




Mass (MeV/c2) Width (MeV/c^) 


a^(1320) 


E852 (PWA+MDF) 
E852 (MDPWA) 


1317 ±1 ±2 127 ±2 ±2 
1313 ±1 119 ±2 


7rf(1400) 


E852 (PWA+MDF) 
E852 (MDPWA) 


1370 ±16 tlo 385 ±40 +^(^5 
1369 ±14 517 ±40 
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Abstract 

Details of the analysis of the r/vr" system studied in the reaction 7r~p -^ r]TT~p 
at 18 GeV/c are given. Separate analyses for the 27 and Tr+vr^vr'^ decay 
modes of the r/ are presented. An amplitude analysis of the data indi- 
cates the presence of interference between the a^(1320) and a J = 1 ^ 
wave between 1.2 and 1.6 GeV/c^. The phase difference between these 
waves shows phase motion not attributable solely to the a^(1320). The 
data can be fitted by interference between the a^(1320) and an exotic 
1-+ resonance with M = (1370 ±16 +^°) MeV/c^ and T = (385 ±40 
_IQ^) MeV/c^. Our results are compared with those of other experiments. 



Typeset using REVT^ 



I. INTRODUCTION 

In a previous publication [|1|, evidence was presented for an exotic meson produced in 
the reaction 

'K~p — >■ ri'K~p (1) 

at 18 GeV/c, with the decay mode rj — > 77. The purpose of this paper is to provide details of 
that analysis, to discuss additional analyses of those data, and to give a detailed comparison 
of our results with those of other experiments. We also compare those results with data from 
our experiment on Reaction (|I|) but with the r] -^ tt+tt^tt^ decay. 

The Tj-K system is particularly interesting in searching for exotic (or non-gg) mesons 
because the system has spin (J), parity (P), and charge-conjugation (C) in the sequence 
J^*" = 0^^, 1^^, 2++, 3 ^ ... for £ = 0, 1, 2, 3, ... . (Here d is the orbital angular momentum of 
the rjTT system.) Hence a resonance with an rjir decay mode with odd l is manifestly exotic. p~| 
Having isospin 1=1, such a resonance could not be a glueball {2g,3g, . . .), but it could be a 
hybrid (qqg) or a multiquark (qqqq) state. 

A. Models 



Properties of hybrids and multiquark mesons have been discussed in the framework of 
various models p-p^. Calculations based upon the MIT bag model predict [§|-§[ that an 
1 = 1, 1 ^ hybrid {qqg) will have a mass near 1.4 GeV/c^. On the other hand, the flux-tube 
model 0,§] predicts the mass of the lowest-lying hybrid state to be around 1.8 GeV/c^. 
Characteristics of bag-model 5- wave multiquark states (which would have J^ = 0^, 1^ or 
2+) have been discussed but those for a 1^ state have not. QCD sum-rule predictions ]10 



'^ A qq meson with orbital angular momentum L and total spin S must have P — (—1)^^^ and the 
neutral member of its isospin multiplet must have C ~ (— 1)^+ . A resonance with quantum number 
in the sequence J^*^ = 0^^,0^ ,1 ^,2^ ,3 ^ ■■■ does not satisfy these conditions and must be 
exotic. 



vary widely between 1.0 GeV/c^ and 2.5 GeV/c^. Recently, an analysis of the multiquark 
hybrids has been carried out, based on the diquark cluster model |]TT|; this model predicts 
a lowest-lying isovector J-^*^ = 1 ^ state at 1.39 GeV/c^ with a very narrow width (~ 8 
MeV/c^). Finally, recent lattice-gauge calculations yield mass estimates for a 1""*" hybrid in 



the range from 1.7 to 2.1 GeV/c [|12 



B. Previous Experiments 

Several experiments, prior to the publication of the E852 results [0], had studied the 
rjiT final state, and observed an enhancement in the P wave around 1400 MeV/c^ |ir3|-|T6 
However, they reached conflicting conclusions. 



The 1988 GAMS experiment |T3[ at CERN {tt~p at 100 GeV/c) claimed to find a nar- 
row enhancement in the unnatural parity exchange Pq wave, but found the natural parity 
exchange P+ wave to be "structureless" . The method of analysis and the conclusions were 
seriously disputed by some of the same authors later [|l^] . 



The 1993 VES experiment |]T4[ at Serpukhov (n^N at 37 GeV/c on a beryllium target) 
found enhancement in the natural parity exchange P+ wave and concluded that "the P+ wave 
is small but statistically significant and contains a broad bump" . They made no attempt to 
identify the "bump" with a resonance. 

The 1993 KEK experiment [|1^ (7r~p at 6.3 GeV/c) claimed that "A clear enhancement 
of the P+ wave was observed around 1.3 GeV/c^" but noted that "The phase of the P+ wave 
relative to the D^ wave shows no distinct variation with mass in the analysis region" . They 
therefore made no attempt to offer a resonance hypothesis. 

The 1994 Crystal Barrel experiment [ITH] on pp annihilation at rest concluded that their 



"TT^vr'^r] data may at most accomodate a small amount of featureless ttt] P-wave" 

The first claim for a 1~"^ exotic resonance in the rjir' channel was made by our experiment 
at BNL |jl|. In this paper we present details of the measurements and method of analysis 
used in that earlier letter publication. We note that since the publication of our letter, an 



independent confirmation of our results has come from a new measurement by the Crystal 
Barrel collaboration ||18[|. 



II. EXPERIMENTAL DETAILS 



A. E852 Apparatus 



Our data sample was collected in the first data run of E852 at the Alternating Gradi- 
ent Synchrotron (AGS) at Brookhaven National Laboratory (BNL) with the Multi-Particle 
Spectrometer (MPS) [|l^ augmented by additional detectors. A diagram of the experimen- 
tal apparatus is shown in Fig. |l]. A Cerenkov tagged vr" beam of momentum 18 GeV/c 
was incident on a one- foot long liquid hydrogen target at the center of the MPS magnet. 
The target was surrounded by a four-layer cylindrical drift chamber (TCYL) ||2^ used to 
trigger on the proton recoil of Reaction (|l]), and a 198-element cylindrical thallium- doped 
cesium iodide array (Csl) |^ capable of rejecting events with wide-angle photons. The 
downstream half of the magnet was equipped with six drift chamber modules (DCl-6) | p2|| , 
each consisting of seven planes, used for charged-particle tracking. Interspersed among these 
were: three proportional wire chambers (TPXl-3) to allow triggering on the multiplicity 
of forward tracks; a window-frame lead scintillator photon veto counter (DEA) to ensure 
photon hermeticity; a scintillation counter (CPVB) to veto forward charged tracks for neu- 
tral triggers; and a window-frame scintillation counter (CPVC) to identify charged particles 
entering the DEA. Beyond the magnet were: a newly-built drift chamber (TDX4) consisting 
of two x-planes; two scintillation counters (BV and EV) to veto non-interacting beam tracks 
and elastic scatters respectively; and a 3045-element lead glass electromagnetic calorimeter 
(LGD) pH] to detect forward photons. Further details are given elsewhere 



B. Trigger 



The trigger (see Ref. [Q for details) for Reaction (|1|) required a recoil charged particle in 
the TCYL detector and one charged particle traversing each of the first two TPX chambers. 
In addition, an electronic algorithm [^ coupling energy and position information in the 
LGD calorimeter (an "effective-mass" trigger) was utilized for the purpose of enhancing the 
fraction of ?7's relative to vr'^'s in the sample. A total of 47 million triggers of this type were 
recorded. 

C. Event Reconstruction and Selection 

Of the 47 million triggers, 47,235 events were reconstructed which were consistent with 
Reaction (|l|). These were selected by requiring: 

• topological and trigger cuts including requirements for: 

1. two photons reconstructed in the LGD; 

2. one forward track reconstructed in DCl-6; 

3. one recoil track reconstructed in TCYL; 

4. a common vertex, in a target fiducial volume, reconstructed from the charged 
tracks and the beam track; 

5. no energy deposited in the DEA detector or outside the fiducial volume of the 
LGD; 

6. the energy deposited in the Csl arrary being less than 160 MeV; 

• that the effective mass of the two photons be consistent with the rj effective mass with 
a confidence level greater than 10~^; 

• that all data come from runs which had proper functioning of the trigger processor; 



that the photons hit the LGD within a fiducial volume which excluded a 4.0-cm region 
(one block width) around the periphery of the LGD as well as a 4.0-cm wide region 
surrounding the beam hole; 

that the distance between a photon and a charged track hitting the LGD exceed 20 
cm; 

-LO < (missing mass)^ < 2.5 (GeV/c^)^; 



• a SQUAW [^ kinematic fit (requiring energy and momentum conservation) to Reac- 
tion (|1]) with a confidence level > 10%; 

• that the difference in angle A0 between the fitted proton direction and the measured 
track in TCYL be less than 8°; 

• the exclusion of events for which the vr" went through an insensitive region of TPX2; 

• the exclusion of events for which the tt~ went through a small region surrounding the 
EV/BV veto counters. (Events which had a vr" in this region were sometimes vetoed, 
probably due to Cerenkov radiation in the EV or BV light pipes). 

Shown in Table | is the effect on the data sample for each of these cuts. The last two cuts 
listed in the table were additional cuts made on the data to carry out the PWA. 

D. Experimental Acceptance 

The experimental acceptance is determined by a Monte Carlo method. Events are gen- 
erated using SAGE |2^ with peripheral production (of the form dN/d\t\ = Aexp —b\t\ with 
b = 4.0(GeV/c)~^) and with isotropic angular distributions in the Gottfried- Jackson (GJ) 
frame. (The GJ frame is a rest frame of the rjir^ system in which the z-axis is in the direc- 
tion of the beam momentum, and the y-axis is in the direction of the vector cross-product 
of the target and recoil momenta.) After adding detector simulation using GEANT |^, the 



Monte Carlo event sample is subjected to the same event-selection cuts and run through the 



same analysis as the data. A second method (called SAGEN) which did not use GEANT 
was employed as well. This method also used SAGE as the event generator, but instead of 
using GEANT, the acceptance was determined using a full detector simulation but without 
such effects as multiple scattering, pair production and secondary interactions. This second 
method allowed acceptances to be calculated much more quickly. The only differences noted 
in the amplitude analysis results (discussed below) between the two methods was in the 
number of events, since GEANT takes account of pair production and secondary interac- 
tions, whereas SAGEN does not. The average acceptances are the ratios of the generated 
events to the accepted events and are shown in Figs. ^-^ using the SAGEN method. 

The average acceptance as a function of r^vr" effective mass is shown in Fig. ^ The 
average acceptance decreases by about a factor of two over the effective mass region from 1.0 
to 2.0 GeV/c^. Average acceptances are calculated for peripheral production and isotropic 
decay as described above. 

Shown in Figs. |] and |^ is the acceptance as a function of cosO and of for various ranges 
of the riTT~ effective mass. Here 6 and (f) are the polar and azimuthal angles measured in the 
GJ frame. The polar angle is the angle between the beam direction and the rj direction in 
this frame. The inefficiency in the backward direction corresponds to slow r^'s and fast tt^'s 
in the lab. The slow ?7's lead to low energy 7's which are often produced at wide angles and 
thus miss the LGD. In some cases, fast tt^'s cause the event to be vetoed if they strike the 
EV or BV scintillation counters, leading to further inefficiency in the backward direction. 
The inefficiency in the forward direction is due to an inefficiency in detecting slow, wide- 
angle pions which can scatter in the Csl detector. The acceptance in is relatively uniform. 
There is a correlation between the energy of the 77 and at finite momentum transfer and 
this leads to the observed shape. 

Finally, shown in Fig.|^ is the average acceptance (GEANT-based) as a function of |t'| = 
\t — tminl; whcre t is the the four-momentum-transfer between the initial- and final-state 
protons and tmin is the minimum value for this quantity for a given tjit' effective mass. 
The dramatic decrease in acceptance below about |t'| = 0.08(GeV/c)^ is due to a trigger 
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requirement. In particular, since we require the presence of a recoil proton in TCYL, the 
trigger cannot be satisfied if the proton stops in the hydrogen target. 

E. Background Studies 

Shown in Fig. |^ is the 27 effective mass distribution for events in the a^(1320) mass region. 
The data sample used for this distribution consisted of a subset of events satisfying the cuts 
listed in Table | but without SQUAW confidence level cuts. The central cross-hatched region 
in Fig.p shows the events which remain after the SQUAW-based kinematic-fitting cuts. The 
distribution has a ~ .03 GeV/c^. Both this distribution and the missing-mass-squared 
distribution discussed below are consistent with that expected from Monte Carlo studies 
when the energy resolution of the LGD for a photon of energy E is taken to be of the form 
a / E = a + hj \fE with a = 0.032 and b = .096 (GeV)^/^. (This was the resolution function 
used for the LGD in the kinematic fitting.) Two methods have been used to study the 
background in our sample. Method 1 used the shaded sidebands of Fig.H and allows us to 
study the non-r^ background in the data. 

The missing-mass- squared distribution for the data sample before kinematic fitting is 
shown in Fig.^ The dashed histogram shows the events which remain after kinematic 
fitting. The distribution for good events for Reaction (|lD should peak at the square of the 
proton mass or at a value of 0.88 (GeV/c^)^ . 

A scatterplot of the missing-mass squared versus the 27 effective mass is shown in Fig. ^ 
Background studies using our Method 2 take as the background estimator a region surround- 
ing the central signal region seen here instead of using the sidebands of Fig. ^ In this way 
we take into account background events of both the non-r^ type and of the type which does 
have an t] present but is not exclusively Reaction || such as events with an extra vr*^. (The 
background is estimated using the region included within the outer elliptical area of Fig.|| 
but not within the middle elliptical region.) 

Shown in Fig.P is the effective- mass distribution of the background, estimated using 



Method 1. In this figure are shown the effective-mass distribution for each side-band region 
as well as the summed distribution for the background regions. Because the background 
regions have different thresholds, one higher than the signal region and one lower than the 
signal region, the histograms are shifted by an appropriate amount (so that their thresholds 
match that of the signal region) before summing. 

In Fig. |10| is shown the polar angular distribution of the background events from Method 
1 in the a^(1320) effective-mass region. Distributions are shown separately for the low-mass 
and the high-mass sidebands of the 77. This high-mass sideband distribution is somewhat 
peaked in the backwards direction with a tendency for the distribution to have an excess 
below the region cos^ < —0.5. We note that this is in the opposite direction from the asym- 
metry in the data (see below) and therefore cannot be the cause of the observed asymmetry. 



Of course the intensity of the background is quite small as seen in Fig. 11 below and therefore 



could not lead to significant changes in the angular distributions in the data in any case. 



III. GENERAL FEATURES OF THE DATA 



The a^(1320) is the dominant feature of the r/vr" effective- mass spectrum shown in Fig. |TT. 
The background, which is shown shaded in the figure, is estimated from Method 2 above, 
and is approximately 7% at 1.2 GeV/c^, and only 1% at 1.3 GeV/c^. 

The acceptance-corrected distribution of |t'| = |t — tmini is shown as the sohd points 
in Fig.O for |t'| > 0.08(GeV/c)^. (Our acceptance is quite low below 0.08 (GeV/c)^ as 



discussed in Section [11 D| .) Since the data are dominated by 02^(1320) production, we show 
as the solid curve (1) the prediction of a Regge Pole model for the differential cross section 
for the reaction 7r~p -^ a^p at 18 GeV/c. (Note that the ordinate values are given by the 
theory and the data are normalized to the theory so we are comparing only the shape of 
the data with the theoretical prediction.) The model includes contributions from p and /2 
Regge trajectories with parameters from a fit by Sacharidis [^ and also includes a small 



at t' = 0.15(GeV/c)^) contribution from a uniform {f independent) background. For 



comparison, results of another experiment (open circles) ||29[ which studied the reaction 
TT^p — > a2P, a2 —* K^K^ at 22.4 GeV/c and the Sacharidis fit (curve 2) are shown. (Note 
that the values shown take into account the a2 decay branching fractions.) We conclude that 
the shape of our t' distribution is consistent with previous experiments and with natural- 
parity exchange production in Regge-pole phenomenology M 



Acceptance-corrected distributions of cos 9 are shown in Fig. |I3| for various ranges of 7]7i~ 
effective mass. For illustration purposes, the acceptance correction is calculated here for 
isotropic decay of the t]it~ system. The acceptance correction used in the amplitude analysis 
discussed below is based upon the observed decay angular distribution. The presence of a 
significant forward-backward asymmetry in the cos 6' distribution is obvious. 

The forward-backward asymmetry in cos 6 is plotted as a function of rjir' effective mass 



in Fig. |T^ Here, the asymmetry is defined as {F — B)/(F + B) where F{B) is the number 
of events in the mass bin with the rj decaying forward (backward) in the G J frame. For this 
figure, the asymmetry was calculated for events in the region with | cos^| < 0.8 to avoid any 
possibility of having results distorted by the extreme forward and backward regions which 
have low acceptance.!^ The asymmetry is large, statistically significant and mass dependent. 
Since the presence of only even values of i would yield a symmetric distribution in cos^, 
the observed asymmetry requires that odd-£ partial waves be present and that they interfere 
with even-£ partial waves to describe the data. Note that the decrease in asymmetry in the 
1.4 GeV/c^ region can be (and will be shown to be) caused by the phase difference between 
the even-i and odd-i waves approaching tt/2 rad. 



The azimuthal angular distribution as a function of r^vr effective mass is shown in Fig. 15. 



The observed structure has a clear sin component which indicates the presence of m = 1 



natural-parity-exchange waves in the production process. (See the discussion in Section |V 
below.) 



Shown in Fig.|16| are the n p and rjp effective-mass distributions for the data sample. 



'' The asymmetry function was plotted for various ranges of the decay angle and the presence of a strong 
asymmetry was noted in all cases. 



It is important to note that the absence of baryon isobar production is required for the 
assumptions of our PWA to be vahd. There is at most a very small amount of isobar 
production in the region M{tt^p) < 2.0 GeV/c^ in Fig.|T^ and none in Fig.|TB|b. The 



amplitude analysis described in Section |^ was checked to insure that isobar production did 
not effect our results. This was done by redoing that analysis after requiring M{7[^p) > 
2.0 GeV/c^ . The resulting intensities and phases did not change (other than an overall 
magnitude change due to the loss of events) in most cases by more than one standard 
deviation and in no case by more than 1.5 standard deviations. 

IV. PARTIAL- WAVE ANALYSIS 
A. Procedure 



A partial-wave analysis (PWA) |[31|-|33| based on the extended maximum likelihood 



method has been used to study the spin-parity structure of the r^vr^ system. We give in 
Appendices A and B some mathematical details regarding the techniques used in the partial- 
wave analysis. The formalism adopted in this analysis is somewhat different from those used 
by previous investigators. Although complete details used in the formalism are given in a 
recent publication by S. U. Chung [^, a portion of that work is reproduced in Appendices 
A and B in order to make this paper as complete and self-contained as possible. 

In Appendix A, a brief description of the formalism is given as are the relationships 
between the partial wave amplitudes (assuming £ < 2) and the moments of the angular 
distribution. The technique of the extended maximum likelihood analysis is covered in Ap- 
pendix B, where the interplay of the experimental moments and the acceptance is described. 
(The experimental acceptance is incorporated into the PWA by using the accepted Monte 
Carlo events described above to calculate normahzation integrals - see ref. pT|). 

The partial waves are parameterized in terms of the quantum numbers J^'" as well as m, 



the absolute value of the angular momentum projection, and the reflectivity e [^. In our 



naming convention, a letter indicates the angular momentum of the partial wave in standard 
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spectroscopic notation, while a subscript of means m = 0, e = —1, and a subscript of +(— ) 
means m = 1, e = +1(— 1). Thus, Sq denotes the partial wave having J^'-^nf = 0^+0^, while 
P signifies l""^!", D^ means 2++1+, and so on. 

We consider partial waves with m < 1 in our analysis. This assumption is true in the 
limit of — t = 0, since the nucleon helicities give rise to the states with m = or m = ±1 
only. But this assumption can be dealt with — experimentally — since the moments H{LM) 
with M = 3 or M = 4 can be checked, to see how important the states \im) are in the data 
with \m\ > 2. This has been done with our data. The moments H{33), H{A3) and if(44) 
are all small in the a^(1320) region, and a fit including |22) (shown in Section [IV Q below) 
contains only a very small amount of this wave and is very broad. 

We also assume that the production spin-density matrix has rank one. This assumption 
is discussed in Appendix C. 

Goodness-of-fit is determined by calculation of a x^ from comparison of the experimental 
moments with those predicted by the results of the PWA fit. A systematic study has been 
performed to determine the effect on goodness-of-fit of adding and subtracting partial waves 
of J < 2 and M < 1. We find that although no significant structure is seen in the waves 
of negative reflectivity (see below), their presence in the PWA fit results in a significant 
improvement in goodness-of-fit compared to a fit which includes only the dominant positive- 
reflectivity partial waves. We have also performed fits including partial waves with J = 
3 and with J = 4. Contributions from these partial waves are found to be within one 
standard deviation of zero for most mass bins with M{ri7T^) < 1.8 GeV/c^ and in all cases 
within two standard deviations of zero. Thus, PWA fits shown or referred to in this paper 
include all partial waves with J < 2 and m < 1 (i.e. Sq, Pq, P-, Doi -D-, -P+, and -D+). A 
non-interfering, isotropic background term of fixed magnitude determined as described by 



Method 2 in Section II D is used. 
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B. Results 

The results of the PWA fit of 38,200 events in the range 0.98 < M{r]Ti-) < 1.82 GeV/c^ 
and 0.10 < |t| < 0.95 (GeV/c)^ are shown in Figs.O and 0. In Fig.O the acceptance- 



corrected numbers of events predicted by the PWA fit for the D+ and P+ intensities and the 
phase difference between these amphtudes, A$, are shown as a function of M{ri7T^). (The 
smooth curves shown in this figure are discussed below in Section [VB| .) There are eight 
ambiguous solutions in the fit [l32| , ^|36| . These solutions are mathematically discrete but 
with equal likelihoods - that is, they correspond to exactly the same angular distributions. 
We show the range of fitted values for these ambiguous solutions in the vertical rectangular 
bar at each mass bin, and the maximum extent of their errors is shown as the error bar. 
These rectangular bars are quite small and thus not apparent for the Z)+ intensity, but they 
are quite clear for the P^-intensity and the phase-difference distributions. 



The a^(1320) is clearly observed in the D^ partial wave (Fig.|T^). A broad peak is seen 
in the P+ wave at about 1.4 GeV/c^ (Fig. |l^b). The phase difference A$ increases through 
the a^(1320) region, and then decreases above about 1.5 GeV/c^ (Fig.[T7|c). This phase 
behavior will allow us to study the nature of the P+ wave. (We note that there is a sign 
ambiguity in the phase difference and thus only the magnitude of A$ is actually measured.) 

Shown in Fig.0 are the fitted intensities for waves which are produced by negative- 
refiectivity (or unnatural-parity) exchange. The predicted numbers of events for these waves 
are generally small and are all consistent with zero above about 1.3 GeV/c^. Although there 
is some non-zero contribution from the D_ and (especially) the Sq waves below this region, 
the uncertainties and ambiguity ranges associated with these waves make it impossible to 
do a definitive study of them to determine their nature. 

The forward-backward asymmetry noted earlier is due to interference in the natural- 
parity exchange sector rather than to the unnatural-parity exchange waves. This is illustrated 



in Fig.lT^ which shows the predicted asymmetry separately for the natural and unnatural 



parity exchange waves. It is clear that the asymmetry due to the unnatural-parity waves is 
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about an order of magnitude less than that due to the natural-parity waves. Also shown in 
Fig. ITUI is the comparison of the asymmetry present in the data with that predicted by the 



fit. The fit clearly does an excellent job in representing the data points. 

The unnormalized spherical harmonic moments H{LM) and their prediction from the 
PWA fit as a function of mass are shown in Fig. |0[ Here H{LM) = J2iLi Yl'^{Oi, 0j) (A^ being 
the number of events in a given bin of M{r]TT^)). The relationships between the moments 
and the amplitudes are given in Appendix A, Eqs. |A27| . An examination of the H{30), 
if (32), H{4:0) and H{A2) moments along with a comparison with Eqs. [A27| shows that the 



D+ amplitude dominates and demonstrates clearly that the P+ partial wave is required for 
the PWA fit to describe the experimental moments. These moments cannot be described 
solely by the combination of the D-|. partial wave and experimental acceptance. 

The change in — log(Likelihood) as a function of the number of events for the P4. partial 
wave for the 1.30 < M{ri7i^) < 1.34 GeV/c^ bin is shown for all the ambiguous solutions in 



Fig. ^ The curves were obtained by fixing the P+ intensity at various values and maximizing 
the liklihood function varying all of the other parameters. All eight ambiguous solutions were 
found for each value of the P+ intensity. The number of predicted P+ events at the maximum 
of the liklihood ranges from 330 events to 530 events for the eight solutions with typical errors 
of 280 events. (A change in — log(Likelihood) of 0.5 corresponds to one standard deviation.) 
For all solutions, the liklihood function gets so bad below 100 events that the P^ wave is 
clearly required to fit the data. Thus the observed variation in — log(Likelihood) further 
demonstrates that the P+ partial wave is required to describe our data. 

C. Systematic Studies 

PWA fits were performed for two different t ranges containing approximately equal num- 
bers of events. One bin spanned the range 0.10 < |t| < 0.25 (GeV/c)^, and the other was for 
0.25 < |t| < 0.95 (GeV/c)^. Both bins yielded comparable structures in the P^ wave, and 
the 02^(1320) was the dominant feature of the D+ wave for both bins. The relative P+-D+ 
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phase behavior for each bin was similar to the resuhs for the integrated fit shown in Fig. |T7|c. 

A PWA fit has been carried out excluding those events with \cos9g.j\ > 0.8 (the region 
in which experimental acceptance is poorest). Neither the P+-wave intensity nor its phase 
variation relative to the D+ wave change by more than one standard deviation in any mass 
bin. 

As mentioned above, a PWA fit has been carried out including the natural parity exchange 
m = 2 amplitude (labelled D2+) in the fit. The results of this fit are shown in Fig.p2|. 
Comparing this fit with the fit shown in Fig. |1^ it is clear that the magnitude and phase 
behavior of the P+ wave is quite unaffected by inclusion of the m = 2 amplitude. 

To test for A and A^* contamination, a fit has been done in which events with M(7r^j9) < 
2.0 GeV/c^ are excluded. As discussed earlier, the resulting intensities and phases did not 
change in most cases by more than one standard deviation and in no case by more than 1.5 
standard deviations. 

Fits were also carried out on Monte Carlo events generated with a pure D^ wave to 
determine whether structure in the P+ wave could be artificially induced by acceptance 
effects, resolution, or statistical fluctuations. Shown in Fig.^ are the results of such a fit. 
We do find that a P+ wave can be induced by such effects. This 'leakage' leads to a P+ wave 
that: (1) mimics the generated D^ intensity (and in our case would therefore have the shape 
of the 02(1320)); and (2) has a phase difference A$ that is independent of mass. Neither 
property is present in our study and we conclude that the P+ structure which we observe is 
not due to 'leakage'. 

Fits have been performed allowing i = 3 and i = A waves. We find that these waves are 
negligible in the region below M{riTT^) of about 1.7-1.8 GeV/c^, their intensities being less 
than one standard deviation from zero in almost all bins. (The largest number of events in 
any bin for the F+ wave was 34 ± 22 events and for the G+ wave was 140 ± 100 events.) 

The data have been fit using different parametrizations of the background. The back- 
ground has been set at fixed values determined from the two different background estimates 
discussed previously. In another fit, the background has been set to zero. And finally, a fit 
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was performed allowing the background level to be a free parameter. Although the nega- 
tive reflectivity waves do change somewhat for different treatments of the backgroundp^, the 
results for the D+ and P+ waves and their relative phase do not change by more than one 
standard deviation in the entire region between 1.2 and 2.0 GeV/c^ except for a few isolated 
points which vary up to 1.5 standard deviations. 

D. Comparison with Previous Experiments 

These results for the P+ and D+ intensities and their phase difference are quite consistent 



with the VES results |jT^ as can be seen in Fig. ^ In particular, the behavior of the shape 
of the phase difference is virtually identical to that reported by VES.Q This is particularly 
noteworthy since, as will be seen below, it is this phase difference which allows us to draw 
conclusions regarding the nature of the P+ wave. 

Our results are compared with those of the KEK experiment iT^ in Fig.pBI. In this case. 



it is clear that the two results differ. The KEK results have a P-wave intensity which is 
narrower and a P-D phase difference which, within errors, is consistent with being constant 
as a function of M{ri7T~). 

E. Comparison ^vith rj -^ 7r+7r~7r*' Data Sample 

A second data set in another topological class (with two additional charged particles in 
the final state) has been used to study Reaction (|1]) with the decay mode rj -^ n^n^iT^. 
Besides having three forward charged particles instead of one, these events have a vr" instead 
of an 1] to be detected by the LGD. Since the tt^ is one of three pions in the r] decay, its 
energy will be significantly less than the t] in the topology with only one charged track. Thus 



'^ Since the background and the Sq wave are both isotropic, the fitting program cannot distinguish 
between them. 

'^ The magnitude of the phase difference is shifted by about 20° relative to that of VES. A produc- 
tion phase shift would not be unexpected because of the differing energies and targets in the two 
experiments. 
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it is clear that the r] -^ ti^ti^ii^ data sample will have significantly different acceptance and 
systeniatics when compared to the 77 -^ 27 sample. 

Shown in Fig.|2^ is the tt+tt^tt^ effective mass distribution from this data set. There is a 
clear r) peak as well as a strong peak in the u region. After kinematic fitting, a sample was 
obtained of 2,235 events which were consistent with Reaction |^ with rj —>■ tt+tt^tt'^. Fig.|27| 
shows the effective mass distribution for this sample of events. As expected, the a^(1320) 
dominates this spectrum. 

Although the data sample is very small for this topology, we have carried out an amplitude 
analysis in order to compare with the primary rj ^ 2'-/ analysis. Results of the analysis are 



shown in Fig. ^ where we compare the shapes of the P+ intensities for the two data sets 
as well as the P+ — D^ phase differences. Despite the rather large statistical uncertainties, 
there is excellent agreement between these distributions. 

V. MASS-DEPENDENT FIT 

In an attempt to understand the nature of the P+ wave observed in our experiment, 
we have carried out a mass- dependent fit to the results of the mass- independent amplitude 
analysis. The fit has been carried out in the r^vr" mass range from 1.1 to 1.6 GeV/c^. In this 
fit, we have assumed that the Z}+-wave and the P+-wave decay amplitudes are resonant and 
have used relativistic Breit-Wigner forms for these amplitudes. 

A. Procedure 

We shall use a shorthand notation w to stand for the rjiT^ mass, i.e. w = Mlrjn"). 
Representing the mass-dependent amplitudes for D+ and P+ as Ve{w) for £ = 2 and 1, we 
may write 

Ve{w) = e'^' Mw) B,{q) [a, + h{w - «;°) + ce{w - w'^.^f^ (2) 

where q is the r/vr" breakup momentum at mass w. Here (pi is the production phase (mass 
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independently) associated with a wave i. The quantities Aii^w) and -B^(g) are the standard 
relativistic Breit-Wigner form and the barrier factor, respectively, and are given below. The 
square-root factor has been introduced primarily to take into account possible deviations 
from the standard Breit-Wigner form, at values of w away from the resonance mass (denoted 
by w^). The overall normalization of a wave is governed by ae, while the constants b^ and q 
allow for deviations in the mass spectra from the Breit-Wigner form. The constants a^, bi 
and Q are all real, so that the square-root factor does not affect the rapidly varying phase 
implied by the standard Breit-Wigner form.p~| 



The barrier functions |37], which are real, are given by: 



Bo{q) 
Bi{q) 

B2{q) 
B-M 
BM 
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1/2 


z + l\ 


Z^ 


1/2 







(3) 



z(z-15)2 + 9(2^-5) 

z^ 



1/2 



[z"^ - Ahz + 105)2 + 2hzi2z - 21)2 



1/2 



where z = {q/q^Y and q^ = 0.1973 GeV/c corresponding to 1 fermi. Note that B£{q) ^ 1 
as g ^ oo. 

The relativistic Breit-Wigner functions can be written 



A Aw) 






TAw) 



'''^("') sin S Aw) 



(4) 



where F^ is the nominal width (mass independent) and r^(w) is the mass-dependent width 
given by 

'BAq)'' 



TAw) = r° 



We 



w 



.4, 



BAqI) 



(5) 



'^ We have tried a linear dependence in mass for the production phase; the fits did not require it. 

^ For a Breit-Wigner form with a constant width, the phase rises 90 degrees over one full width centered 
at the resonance mass. 
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where g° is the breakup momentum evaluated at w = w^. The mass- dependent phase shift 
Si{w) is given by 



cot 6e{w) 



TAw) 



/ w 



(6) 



or 



-^^^-)-te)(7 



BM) 



-\2 



and the overall phase for the £-wave amphtude is 



^E = 4>i + 6i{w). 






(7) 



We are dealing with two waves, P+ and -D+, and can only measure 
phase difference being measured experimentally, corresponds to 

A$ = $2 - $1 = (/) + 52iw) - 5i{w) 



()2 — 01. Thus the 



(9) 



Finally, the experimental mass distribution for each wave £ is given by 

^ = \V,{w)\'pq 



(10) 



where pq is the phase-space factor for which p is the breakup momentum of the r/vr^ system 
(or of the final-state proton) in the overall center-of-mass frame in Reaction (|T]). Since the 
problem here is for a given ^/s, all other relevant factors, including that of the beam flux, 
have been absorbed into the amplitude itself, i.e. the constants ae, be and q. 

The input quantities to the fit included, in each mass bin: the P+-wave intensity; the 
D_i_-wave intensity; and the phase difference A$ (the relevant formulas are given in (^ and 
(p!0|)). Each of these quantities was taken with its error (including correlations) from the 
result of the amplitude analysis. One can view this fit as a test of the hypothesis that the 
correlation between the fitted P-wave intensity and its phase (as a function of mass) can be 
fit with a resonant Breit-Wigner amplitude. 



We find that the fit does not improve significantly when the P+ wave is modified from 
the Breit-Wigner form, and hence set bi and Ci = for the final fit. We also note that the 
magnitudes of the quantities 62 and C2 in the final fit correspond to a small deviation of the 
L'+'Wave intensity of the order of 1%. 

B. Results 

Results of the fit are shown as the smooth curves in Fig. p!7|a, b, and c. The mass and 
width of the J^^ = 2++ state (Fig. 0a) are (1317 ±1 ±2) MeV/c^ and (127 ±2 ±2) MeV/c^ 
respectively |3^. (The first error given is statistical and the second is systematic.) The mass 



and width of the J^'^ = 1"+ state as shown in Fig. 0b are (1370 ±16 +3J) MeV/c^ and (385 
±40 J^^pg ) MeV/c^ respectively. Shown in Fig. p!7| d are the Breit-Wigner phase dependences 
for the a^(1320) (line 1) and the P+ waves (line 2); the fitted D+— ?+ production phase 
difference (line 3); and the fitted D+— P+ phase difference (line 4). (Line 4, which is identical 



to the fitted curve shown in Fig. p!7|c, is obtained as line 1 — line 2 ± line 3.) 

The systematic errors have been determined from consideration of the range of solutions 
possible because of the ambiguous solutions in the PWA. Since there are 8 ambiguous solu- 
tions per mass bin and we are fitting over 12 mass bins, it is clearly impossible to try all 8^^ 
possible combinations. Instead, we have fit some 10^ combinations where the values to be 
fitted in each mass bin have been chosen at random from among the 8 ambiguous PWA solu- 
tions. The resulting fits generally clump into a group with reasonable values of x^ /dof{< 2) 
and into a group with poor values. The systematic errors on the mass and width given above 
are taken from the extremes observed for the solutions with reasonable values of x^ /dof . 
The central values quoted above are taken from a fit which uses the average values of the 
input parameters in each bin. 

The fit to the resonance hypothesis has a x^/dof of 1.49. The fact that the production 
phase difference can be fit by a mass-independent constant (of 0.6 rad) is consistent with 
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Regge-pole phenonienologyp~| in the absence of final-state interactions. If one attempts to 
fit the data with a non-resonant (constant phase) P+ wave, and also postulates a Gaussian 
intensity distribution for the ?+ wave, one obtains a very poor fit with a x^/dof of 7.08. 
Finally if one allows a mass- dependent production phase, a x^/dof of 1.55 is obtained for 
the non-resonant hypothesis — but in this case the production phase must have a very 
rapid variation with mass.[^ Furthermore, for this non-resonant hypothesis, one must also 
explain the correlated structure observed in the P+ intensity — a structure which is explained 
naturally by the resonance hypothesis. 

An attempt ||3^ to explain our result as the interference of a non-resonant Deck-type 
background and a resonance at 1.6 GeV/c^ can reproduce this correlation. (Evidence for an 
exotic meson with a mass near 1.6 GeV/c^ has been reported ^^J by our collaboration.) How- 
ever, this explanation is excluded because of the recent observation |18[ by the Crystal Barrel 
collaboration which confirms the presence of this state produced in nucleon-antinucleon an- 
nihilation. The Deck-effect is a mechanism applicable to peripheral production but not to 
annihilation. 

Our fitted parameters for the J^*^ = 1~+ resonance are compared in Table |I| with the 
values reported by the Crystal Barrel experiment |T^. That experiment reports that a 
J^*^ = 1""'" resonance in the rjir channel is required to fit their data in the annihilation 
channel pn -^ 71^71^1]. Their fitted parameters are very consistent with those determined 
from our mass-dependent analysis. 

C. Other Systematic Studies 



^ The signature factor and the residue functions are at most t-dependent (not mass dependent) (see 
ref. 0). 

'* The fit requires a hnear production phase difference with a slope of -4.3 rad/GeV. 
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1. Sensitivity to the D-wave Intensity Distribution Function 

In order to determine the sensitivity of the results of our mass- dependent analysis to 
the exact function being used to fit the D-wave intensity distribution we have redone the 
fit using two other hypotheses. First we have performed a fit in which the mass-dependent 
amplitude is given by Eq.(0), but with 62 = C2 = 0. Second, we have taken 62 = C2 = in 
Eq.(^ and also replaced the Blatt-Weisskopf barrier functions for each wave by the factor 
q^. Although the resulting fits are poorer in quality, we find that the parameters of the fit 
do not change by amounts greater than the systematic uncertainty described above. 

2. Sensitivity to Leakage 

As shown in Section [IV C| , the P+ wave observed in our data is not consistent with 
"leakage" . That is, the analysis shows that the intensity and phase motion of the P+ wave 
do not have the characteristics of the wave which is artificially generated from a pure D^ 
wave due to possible incomplete knowledge of the resolution or detection inefficiency. This 
does not preclude the possibility of some leakage being present in the data and distorting 
the results of the mass dependent analysis (MDA). In this section, we describe a test which 
has been carried out to study the sensitivity of our MDA results to possible residual leakage 
being present in the data. 

The fit which has been carried out is a mass dependent partial wave analysis (MDPWA). 
In such a fit, the PWA is carried out as in Section |^ but instead of carrying it out separately 
for each t]7t~ mass bin, all bins are fit simultaneously and are tied together with a mass- 
dependent function for each partial wave. That is, the extended maximum likelihood function 
of the form given by Eq. [B16| is generalized to include mass dependence: 

n „ 

ln£ oc Y. ln/(fii, Wi) - / dQdw r/(fi, w) I{Q, w). (11) 

i 

The free parameters in the fit include, in addition to the amplitudes of the partial waves, 
the Breit-Wigner masses, widths, and intensities as well as mass-independent production 
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amplitude phases. 

For simplicity, we have taken the Breit-Wigner form of Eq. ^ to describe each of the 
partial waves. A common mass and width were used for the /)+, D_, and Dq partial waves. 
Similarly, the P+, P_, and Pq waves were assumed to be described by a common mass and 
width. The 5* wave was assumed to have its own mass and width. The constants bg and q 
of Eq. 1^ were taken to be zero for the Dq, D^, P_|_, Pq, P_, and 5*0 waves. Each wave was 
allowed to have its own normalization constant and production phase. 

In order to include leakage in the fit, a leakage amplitude P^'^ with the characteristics 
obtained in the leakage study of Section |IVC| was defined. This amplitude was taken to have 



the shape of the D+ amplitude as well as its Breit Wigner phase dependence. Its production 
phase was fixed at (p^^''^ = 80° This amplitude was then combined coherently with the P+ 
signal to give an effective amplitude given by the expression: 

P'f^f\w) = P+{w) + Pf'\w) (12) 

where 

p!1'\w) = af ) e^^'"' A2{w- w^, T^) B,{q) [1 + ht{w - wl) +ht{w - wlfY/\ (13) 



The results of the MDPWA fit are shown as the smooth curves in Fig. ^. Also shown 
as the points with error bars are the results of the mass-independent PWA. It is clear that 
the two analyses give consistent results. Shown in Fig. p9| a are the P'^ intensity (curve 
3) along with the P+ signal intensity (curve 1). The leakage is shown in curve 2 as the 
sum of the leakage intensity and the (signal - leakage) interference term. The fitted leakage 
contribution is equal to TZ^ = 0.018 where we have defined 

n, ^ [^) . (14) 

The results of the fit are given in Table ^ where they are compared with those of the 
combined PWA and mass dependent fit. The results are quite compatible when one takes 
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into account the systematic errors. The biggest difference is in the fitted width of the P+ 
state which is larger for the MDPWA. 

In Fig. pUI are shown the fitted values for the mass (curve 1) and the width (curve 2) of 



the l""*" resonance as well as the change in ln£ as a function of TZi (the leakage fraction). 
We note that the mass and width are very insensitive up to values of TZ^ = 5%, above which 
the fit becomes very unlikely. 

D. Cross Section Estimate 

In order to estimate the cross section for production of the observed P+ state (which we 



now refer to as 7ri(1400)), we have fitted published total cross sections PT|-[5T[] for 02^(1320) 
production to a function of the form a = A{p/pq)~"' + B, where p is the beam momentum and 
Po is set to 1 GeV/c. Experiments with poor a]^(1260)/a2"(1320) separation were excluded 
from the fit. The best fit gave: A = 5099 ± 221 fib; n = 1.88 ± 0.03; and B = 39.2 ± 2.0 fib. 
From this we estimate the total cross section for a^(1320) production at 18.2 GeV/c to be 
61.1 ± 2.2 fib. This is in good agreement with the result of 62.56 ± 2.92 fib measured [^ at 
18.8 GeV/c. 

From the results of our PWA, we find that, in the rj mass range from 1.10 to 1.58 GeV/c^ 
there are 60, 332 ± 2, 060 -0+ events, and there are 3, 321 ± 1, 245 P+ events. Here the error 
for the number of D+ events is statistical and the error for the number of P+ events includes 
uncertainties due to ambiguities. One thus obtains: a^n^p -^ pvrf (1400)) *PP(7rf (1400) -^ 
TfTT-) = 0.49 ±0.19 fib. 

VI. SUMMARY AND CONCLUSIONS 

In this paper, we have discussed the details of the amplitude analysis of data from Reac- 
tion 0. Interference between D-wave and P-wave amplitudes produced with natural parity 
exchange is required in order to explain the data. Using this interference, we have shown 
that the P-wave phase has a rapid variation with mass and that this phase variation coupled 
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with the fitted P-wave intensity distribution is well-fitted by a Breit-Wigner resonance with 
mass and width of (1370 ±16 +^J) MeV/c^ and (385 ±40 ^^^^) MeV/c^ respectively. Since 
a P-wave resonance in the rjiT system has J^'" = 1 *", it is manifestly exotic. The exact 
nature of the observed state awaits further experimentation. 
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APPENDIX A: PARTIAL- WAVE FORMULAS 

In this appendix, the angular distributions are derived for the rjiT^ system produced in 
Reaction ([l|). The distributions are given both in terms of the moments and the amplitudes 
in the reflectivity basis. For a system consisting of S, P and D waves, explicit formulas for 
the moments as functions of the partial waves are also given. 

In the Gottfried- Jackson (GJ) frame, the amplitudes may be expanded in terms of the 
partial waves for the r/vr" system: 

Uk{n) = Y, V^fmfc^m(fi) (Al) 

where Vimk stands for the production amplitude for a state \im) and k represents the spin 
degrees of freedom for the initial and final nucleons {k = 1,2 for spin-nonfiip and spin- flip 
amplitudes). Aim{^) is the decay amplitude given by 



AU^) = ^^ Di*(0, e, 0) = Yrm (A2) 
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where the angles Q = {6, (p) describe the direction of the r] in the GJ frame. It is noted, 
in passing, that the small (i-function implicit in ( |A2D is related to the associated Legendre 
polynomial via 



dmoi^) 



The angular distribution is given by 



^I^^'-f"^'" '^^' 



I{Q)=Y.\Ukm^ (A4) 

k 

It should be emphasized that the nucleon helicities are external entities and the summation 
on k is applied to the absolute square of the amplitudes. A complete study of the rjn' 
system requires four variables: M{r]7c^), —t and the two angles in Q. The distribution ( [A^ ) 
is therefore to be applied to a given bin of M{ri'K^) and of —t. 

The angular distribution may be expanded in terms of the moments H{LM) via 

n^) = E (^^) H{LM)Di* i^, e, 0) (A5) 

where 

/Of 4-]\ ^/^ 

HiLM) = Y: -^^ p'^^,{i'm'LM\im){i'OLOm (A6) 

i'm' 

where p is the spin-density matrix given by 

Pmrri' = E ^ (-rnk^ i' m! k (A7) 



It is seen that the moments H{LM) are measurable quantities since 

H{LM) = Jdn i{n) Dl,Q{(f), e, o) (as) 

The normalization integral is 

H{00) = f dQ I{Q) (A9) 
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The symmetry relations for the if 's are well-known. From the hermiticity of p, one gets 

H*{LM) = (-)*^ H{L - M) (AlO) 

and, from parity conservation in the production process, one finds 

H{LM) = {-Y^ H{L - M) (All) 

These show that the if s are real. 

Parity conservation in the production process can be treated with the reflection operator 
which preserves all the relevant momenta in the S-matrix and act directly on the rest states 
of the particles involved. It is important to remember that the coordinate system is always 
defined with the y-axis along the production normal. In this case the reflection operator is 
simply the parity operator followed by a rotation by vr around the y-axis. 

The eigenstates of this reflection operator are 

\eim) = eim){\im) - e(-)"K - m)} (A12) 

where 

6{m) = —i=, m > 
v2 

= -, m = (A13) 

2' 

= 0, m < 
For positive reflectivity, the m = states are not allowed, i.e. 

\em) =0, if e = + (A14) 

The reflectivity quantum number e has been defined so that it coincides with the naturality 
of the exchanged particle in Reaction (|l|). One can prove this by noting that the meson 
production vertex is in reality a time-reversed process in which a state of arbitrary spin- 
parity decays into a pion (the beam) and a particle of a given naturality (the exchanged 
particle) 
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JVj ^ s^^ + TT (A15) 

where 77's stand for intrinsic parities. The hehcity-couphng amphtude F^ for this decay |^ 



IS 



Al{M) o: Fi Dit,{<P,,e,,Q) (A16) 

where A is the hehcity of the exchanged particle and the subscript p stands for the 'produc- 
tion' variables. M is the z-component of spin J in a J rest frame. From parity conservation 
in the decay, one finds 

F/ = -F\ (A17) 

where one has used the relationships 77^ = {—Y (true for two-pseudoscalar systems) and 
rjs = (— )* (natural-parity exchange). The formula shows that the helicity-coupling amplitude 
F"^ is zero if A is zero. Since angular momentum is conserved, its decay into two spinless 
particles in the final state cannot have M = along the beam direction (the GJ rest system), 
i.e. the D-^-function is zero unless M = A, if ^p = 0p = 0. Finally, one may identify J with 
i and M with m, which proves ( |A14| ). 



The modified D-f unctions in the reflectivity basis are given by 

^^i*o(0, 0, 0) = 9im) [d^*o(0, 0, 0) - 6(-)- ^1^(0, 0, 0)] (MS) 

It is seen that the modified D-functions are real if e = — 1 and imaginary if e = +1: 

(-)D^*(</),^,0) = 29{m)dU0)cosm<f) 

(A19) 

The overall amplitude in the reflectivity basis is now 

'Uk{n)=Y.'Vimk'^im{^) (A20) 

im 
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where 



'AerrM = f-^^Dlti<P,9,0) (A21) 

and the resuhing angular distribution is 

/(f])=5:rf/fc(f])P (A22) 

ek 

It is seen that the sum involves four non-interfering terms for e = ± and k = 1,2. The 
absence of the interfering terms of different reflectivities is a direct consequence of parity 
conservation in the production process. We use the partial wave amplitude notation 

[^]o = ^-V,o, [i]- = ^-^Vn, [i]+ = W\4i (A23) 

where [i] stands for the partial waves S, P, D, F and G for i =0, 1, 2, 3 and 4. 

Consider an example where the maximum £ is 2. One sees that there are in general twelve 
possible non-zero experimental moments: 



iJ(OO), if(lO), H{11), H{20), H{21), H{22) 
H{30), H{31), H{32), H{40), H{41), H{42) 



(A24) 



while the partial waves [i] are, for unnatural-parity exchange, 

^0, Po, P-, ^0, D_ (A25) 

and, for natural-parity exchange, 

P+, D+ (A26) 

One wave in each naturality can be set be real (5*0 and P+, for example), so that there are 
again twelve real parameters (to be determined). It is helpful to write down the moments 
explicitly in terms of the partial waves: 
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H{10) = -^SoPo + -^PoDo + -^{P-D- + P+D+) 

Hi20) = ^SoDo + Ip^ - liP'_ + Pi) + ^-Dl + \{Dl + Dl) 

H{22) = \^-{Pl-Pl) + \J\{Dl-Dl) 

3 (A27) 

^(30) = — =(v^Po^o - P-D. - P+D+) 

H{31) = iy'^(2Po^- + ViP^Do) 

HiAl) = \\JIdoD. 
H{A2) = y^{Dl - Dl 



^l) 



APPENDIX B: MAXIMUM-LIKELIHOOD ANALYSIS 

This appendix is devoted to an exposition of the experimental moments, the acceptance 
moments and the acceptance-corrected (or 'true') moments and the relationships among 
them. Finally, the extended likelihood functions are given as functions of the 'true' and 
acceptance moments. 

One may determine directly the experimental moments (unnormalized) as follows: 

n 

H,{LM) = Y^ Di,,{<P,,,e,,0) (Bl) 

i 

where the sum is over a given number n of experimental events in a mass bin. But this is 
given by, from ( |A8| ) , 
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H^{LM) = Jdn r]{n) i{n) Di,^{<p, e, o) (B2) 

where r]{fl) represents the finite acceptance of the apparatus, and includes software cuts, if 
any. From ([A5[ ), one finds that 

H^{LM) = J2 H{L'M') ^,{LML'M') (B3) 

L'M' 

where 

/or' + i\ r 
^^{LML'M') = ( ^^1 J dfi Vm D^,,,i<p,e,0)D^,',*{<p,e,0) (B4) 

Note that the \l/'s have a simple normalization 

^.,{LML'M') = 5ll'Smm' (B5) 



in the limit ri(il) = 1. The integral ( [B4| ) can be calculated using a sample of 'accepted' MC 
events. Let A^^^ be the number of accepted MC events, out of a total of A^ generated MC 
events. Then, the integral is 

/ 2 /' ' + 1 \ 1 ^"^ 
^^{LML'M') = (^^1 ^E DiU<p,,9,,0)D^^*{cP,,9,,0) (B6) 

Equation ( |B3| ) shows that one can predict the experimentally measurable moments (|Bl|), 
given a set of true moments {H} and the \E''s; this provides one a means of assessing the 
goodness of fit by forming a x^ based on the set {H^}. 

There exists an alternative method of determining \E''s. For the purpose, one expands 
the acceptance function t]{Q) in terms of the ortho normal D-functions, as follows: 

V{n) = E(2^ + ^nLM) Di,t (0, e, 0) (B7) 

LM 

where ^(LM) is given by 

^LM) = ^Jdn 7^{n) D^,o(0, 0, 0) (B8) 
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The complex conjugate is, from the defining formula above, 

C{LM) = {-Y^i{L -M) (B9) 

so that the acceptance function can be made explicitly real 

r^iSl) = Y,(2L + l)r(M)Re [i{LM) Di,*{4>, 9, 0)} (BIO) 

LAI 

where 

t{M) = 2, M > 0, 

= 1, M = 0, (Bll) 

= 0, M < 



One sees that t{M) = 46*2 (M) where e{M) is defined in Eq. i \KT^ 



A set of ^{LM) specifies completely the acceptance in the problem. The normalization 
for the acceptance function has been chosen so that a perfect acceptance is given by 77 (fi) = 1 
and ^(LM) = 6loSmo- The ^(LM)'s can be measured experimentally using the accepted MC 
events 

aLM) = ^^ E ^Mo(0., 0^, 0) (B12) 



Finally, substituting ( p7|) into ( p^ , one finds 

<^^{LML'M')= Y. (21" + 1)C{L"M"){LML"M"\L'M'){L0L"0\L'0) (B13) 

L"M" 

This formula shows an important aspect of the ^{LM) technique of representing acceptance. 
Although (p^) involves a sum in which L and M could be extended to infinity for an arbitrary 
acceptance, there is a cutoff if the set {H} has maxima Lm and M^ [see (p3|) ]. The formula 
above demonstrates that L" < 2Lm and \M"\ < 2Mm- 

In a partial-wave analysis, it is usually best to take a set of the partial waves, [£]o, [C\~ 
and [£] + , as unknown parameters to be determined in an extended maximum-likelihood fit. 
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Since there is an absolute scale in an extended maximum-likelihood fit, one then has the 
predicted numbers of events for all the partial waves, corrected for finite acceptance and 
angular distributions. The partial waves in turn give rise to a set of predicted moments 
{H}. But the H{00) is not 1 but the total predicted number of events from the fit, i.e. one 
should be using the unnormalized moments. It is possible to choose iJ's as unknowns in the 
fit, but the two sets of iJ's should be the same ideally — this affords one an effective way of 
assessing self-consistency between the chosen moments and the partial waves. 

For completeness, a short comment is given about the extended likelihood functions. 
The likelihood function for finding 'n' events in a given bin with a finite acceptance 1]{Q) is 
defined as a product of the probabilities. 



C oc 



n" 



n\ 



n 



ji{n)ri{n) dn 



(B14) 



where the first bracket is the Poisson probability for 'n' events. This is the so-called extended 
likelihood function, in the sense that the Poisson distribution for 'n' itself is included in the 
likelihood function. Note that the expectation value n for n is given by 



noc f I{n)r]{n) dVL 



(B15) 



The likelihood function C can now be written, dropping the factors depending on n alone. 



C oc 



n Ii^^ 



exp 



i{n) 7]{vl) dn 



The 'log' of the likelihood function now has the form, 

n „ 

ln£ oc J2 ln/(fii) - dQ r]{Q) I{Q) 



(B16) 



which can be recast in terms of the ^(LM)'s 

n 

ln£ oc ^ ln/(a) - ^ (2L + 1) H{LM) C{LM) 

i LM 

n 

OC ^ ln/(a) - ^ (2L + 1) t{M) H{LM) Re^(LM) 



(BIT) 



LM 
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H{LM)^s may be used directly as parameters in the fit or may be given as functions of tfie 
partial waves. It is interesting to note that the ^(LM)'s for L > L^ and |M| > M^ are not 
needed in the likelihood fit. Note also that only the real parts of the ^(LM)'s are used in 
the fit. 

It should be borne in mind that a set of the moments {H} may not always be expressed 
in terms of the partial waves. This is clear if one examines the formulas ( |A271 ). Consider, 
for example, an angular distribution in which H{10) is the only non-zero moment. But this 
moment is given by a set of interference terms involving even-odd partial waves. So at least 
one term cannot be zero — for example, the interference term involving S- and P-waves. But 
then neither if (00) nor if (20) can be zero, since both S- and P-waves are non-zero. One 
must conclude then that a x^ based on the set {H^} may not necessarily be zero identically. 

APPENDIX C: RANK OF THE DENSITY MATRIX 

An assumption needed for the partial- wave analysis is that the density matrix has rank 1, 
i.e. the spin amplitudes do not depend on the nucleon helicities. Our justification, so far, has 
been that the fitted partial waves are very reasonable, that these waves can be fitted with 
a very simple mass- dependent formula, that Pomeron-exchange amplitudes are in general 
independent of nucleon helicities, and so on. . . . 

The purpose of this appendix is to point out that, under a simple model for mass depen- 
dence of the partial waves, it is possible to prove that the spin density matrix has rank 1. 
Suppose that one has found a satisfactory fit under a rank-1 assumption. One can then show 
that, even if the problem involves both spin-nonflip and spin-flip at the nucleon vertex — i.e. 
it appears to be a rank-2 problem — the spin density matrix in reality has a rank of 1. Al- 
though this note is based on the results of our r^vr^ analysis, the derivation does not depend 
on the decay channels; the conclusions apply equally well to other decay channels. 

This note relies on some technicalities generally well known, and so they have been 
presented without attribution. The reader may wish to consult a number of preprints and/or 
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papers, which deal with them in some detail pT|j3^,|52|-p4 



1. Partial Waves Produced via Natural-parity Exchange 

Consider the rj7r~ system produced via natural-parity exchange. It consists mainly of just 
two waves -D+ and P+ in the a^(1320) region. Assume these are the only waves. Without 
loss of generality, the decay amplitudes [^ can be considered real, i.e. 



A^{n) = \ —V^dl^i^d) sin0= -J— V3 sin0cos^sin0 
V 47r V 47r 

A^{n) = J^ V2dlf^{9) sin0 = -J^ sin^ sin0 



(CI) 



Since one deals with the partial waves produced only by natural-parity exchange, one can 
drop the subscript '+' from the waves, and the angular distribution is simply given by 



3 \ ^_ „ _ 2 

Air 
3 



oc 



oc 



An 



yhD cos 6 + P sin^ 9 sin^ 

'SIZ^I^ cos'^e + 2V5^{D*P} cos^+|P|2 



(C2) 



sin 9 sin^ 



The integration over the angles can be carried out easily, to obtain 



f I{n)dn(x \D\^ + \P\ 



(C3) 



as expected. 

It is easy to calculate the forward-backward asymmetry A{F, B) = {F — B)/{F + B 
(see section |IlD 

3^5 ,._ \P\\D 



A{F,B) = ^— cos(A$) 



4 



(|P|2 + |D|2) 



(C4) 



where A<l> is the phase difference between the P and D waves. 
The spin density matrix is given by 

I{Q) (X \DA^{Q) + PA^{Q)\^ = ^ Pk,k'AkAl, 

k,k' 



(C5) 
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where {k, k'} = {1, 2} and '1' ('2') corresponds to D (P). From this definition, one sees that 

/ |D|2 DP*\ 

P=[ , J (C6) 

\D*P |P|2 / 

One can work out the eigenvalues of this 2x2 matrix: 

A = {|Dp + |Pp, 0} (C7) 

One of the two allowed eigenvalues is zero, i.e. the rank of this matrix is 1. This is the 
'rank-1' assumption one makes to carry out the partial- wave analysis and is valid for a given 
mass bin. 

Suppose now that the rank is 2, i.e. 

i{n)oc\D,A^{n) + p,A^{n)\^ + \D2A^{n) + P2A^{n)\^ (cs) 

where subscripts 1 and 2 stand for spin-nonfiip and spin-flip amplitudes at the nucleon vertex 
for Reaction ([l|). Comparing (|C2|) and ( |C8|) , one flnds immediately 



|P|2=|p,|2 + |p2|2 (C9) 

^{P* D} = 3fJ{Pi* D^} + ^{P; D^} 

Let w be the effective mass of the r]7T~ system. If the mass dependence is included 
explicitly in the formula, one should write, in the case of rank 1, 

^^^ oc \D{w) A,{Q) + P{w) A,{Q)f pq (CIO) 

where p is the breakup momentum of the "qn^ system in the overall CM system and q is the 
breakup momentum of the rj in the rjir' rest frame. Note that both p and q depend on w. 
Note also that the w dependences of the partial waves D and P are given in the formula. 
Obviously, a similar expression could be written down for the case of rank 2. 
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One is now ready to make the one crucial assumption for a mass-dependent analysis of 
the D and P waves: one assumes that two resonances — in D and P waves, respectively — are 
produced in both spin-nonflip and spin-flip amplitudes. One may then write, for the rank-1 
case, 

D{w) = ae'"e'^^ sin^a 

(Cll) 

P{w) = he'^'' sm5b 

where a, h and the production phase a are all real and independent of the 777?" mass. In 
addition, one can set a > and 6 > without loss of generality. Here 5a and 5}, are the 
phase-shifts corresponding to the resonances and are highly mass dependent. In its generic 
form, the Breit-Wigner formula is given by the usual expression 

cot5 = ^i^ (C12) 

where Wq and Fq are the standard resonance parameters. In this note, the width is considered 
independent of w. Likewise, the barrier factor dependence for D and P is ignored.^ 
The formulas (|C11| ) are generalized to the case of rank 2, as follows: 

Diiw) =aie^"^e^'''^ Sin5a 

P^{w) =6ie*^' sin^fe 

(C13) 

D2H =a2e*"^e*'''" sin^^ 
P^{w)=b2e'^'- sin 6b 

Once again, a^, 6j and a^ are real, a^ > and hi > 0, and independent oi w. One finds, using 



Although simphfied formulas are used in this note, the results given here do not change even when 
correct formulas are used. Note that, to go over to a correct formulation for each wave, one needs to 
substitute the absolute value of the Breit-Wigner formula as follows: 



sin(5(it;) -^ B{q) 



T{w) 



sm5{'w) 



where B{q) is the barrier factor and V{w) is the mass-dependent width. It should be noted that the 
correction factors are all real, by definition. 
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2 2 2 

a = ai + a2 

b-' = bl + hi (C14) 

a6cos(a + 5a- Sb) = cti&i cos(q;i + 6a - Sb) + 0262 cos(a2 + Sa - Sb) 

A plot of cos(a + Sa — 6b) as a function of w is shown in Fig. ^ for three values of a, i.e. 
0°, 45° and 90°. The resonance parameters for a and b of 1.0 and 0.151 are taken from the 
mass dependent fit of Section |V B| as are the resonant masses and widths.^ The normalized 



absolute squares of the Breit-Wigner forms are given in Fig. 32, as is the 'normalized' inter- 



ference term. The same quantities, as they appear in Ref. [|I[, are shown in Fig. ^. This 
figure shows how important the interference term is compared to the P-wave term. Note 
also how rapidly the interference term varies as a function of w in the 02 (1320) region. This 
term, of course, is intimately related to the asymmetry in the Jackson angle and vanishes 



when integrated over the angle, i.e. it does not contribute to the mass spectrum [see (|C^ ) 



and (P3|)]. Fig. |3^ shows the contour plot of the intensity distribution in w vs. cos 6'; note 
the variation of the asymmetry as a function w. 

For the last equation in ( |C14| ) to be true for any mass, the coefficient of cos{6a — Sb) or 



sin(5a — 6b) on the left-hand side must be equal to that on the right-hand side, so that 

ab cos a = ai&i cos ai + 02^2 cos 02 

(C15) 

ab sin a = aibi sin ai + 0262 sin 0^2 

Taking the sum of the squares of the two formulas above and introducing the first two 
equations of (|C14| ), one obtains: 



2ai6ia2&2 cos ai cos 0:2 + 2ai6ia2&2 sin ai sin 02 

--ajbl + ajbj (C16) 

=0^62 (cos^ ai + sin^ ai) + a\b\ (cos^ 0^2 + sin^ 0^2) 



^ The value of a as determined from this fit is 37.46°; for the purpose of illustration, one may consider 
a = 45° close enough. 
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which is recast into 

= (0162 cosai — 0261 cos 0:2)^ + (ai&2 sinai — 0261 sin 02)^ (C'17) 



It is clear that each term must be set to zero, so that 

ai\ fO'2\ 

-— cosai =1— cos 02 

biJ V 62 / 

ai\ . fa2\ . 

-— sm «! = -— sm ^2 
hi J V 02 / 



(C18) 



Placing these back into (|C15| ), one deduces that 



a\ fai\ (a2 

— I cos ct = I 7- ) cos c^i = I 7- ) cos 0:2 

a\ . /ctA . (0.2 



(C19) 



One may take — alternately — the sum of the squares of the two formulas above, or a division 
of the second by the first, and obtain (remembering that the a's and 6's are non-negative 

real quantities), 

a Oi 02 

h h &2 (C20) 

tan a = tan ai = tan 02 

The last equation above demands that ai and ^2 are determined (up to ±7r), but they have 
to satisfy (|C19|) . It is therefore clear that one must set a = ai = 02. Next, one introduces 
two new real variables a; > and y > 0, given by 

(C21) 



X 

a b 

02 62 

a 



with the constraint x"^ + y"^ = 1. 

Now one can prove that the case of rank 2 is reduced to that of rank 1. Indeed, one sees 
immediately that 



and ( |C8| ) becomes identical to ( P2| 
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2. Discussion 

It is shown in this appendix that the problem of two resonances in D^ and P+ in the r/vr" 
system in (|1|) is — effectively — a rank-1 problem. For this to be true, the following conditions 
have to be met: 

(a) There exist two distinct resonances with different masses and/or widths. Note that the 
crucial step, from (|C14|) to (|C15| ), depends on that fact that 5a — 5f, is non-zero and is 



mass dependent. 

(b) There exists a satisfactory rank-1 fit with two resonances in a given mass region, in 
which each amplitude for D+ or P+ has the following general form 

7Wfc(w;,fi) = rfce*°^e^^'=("')/,HAfe(f]) (C23) 

where k = {1,2} and '1' ('2') corresponds to D^ (-P+)- 5k{w) is the Breit-Wigner 
phase and highly mass dependent, while r^ and a^ are mass independent in the fit. 
Of course, one of the two a^'s can be set to zero without loss of generality, so that 
there are three independent parameters, e.g. ri, r2 and ai (these were denoted a, b 
and a, respectively, in the previous section). fk{w) contains the absolute value of the 
Breit-Wigner form, plus any other mass- dependent factors introduced in the model. 
y4fc(n) carries the information about the rotational property of a partial wave k. 

(c) The same two -D+ and P+ resonances are produced in both spin-nonflip and spin-flip 
amplitudes, with the same general form as given above — but with arbitrary r^'s and 
afc's for each spin-nonflip and spin-flip amplitude. It has been shown in this appendix 
that only one set of r^'s and a^'s, i.e. ri, r2 and ai, is required for both spin-nonflip and 
spin-flip amplitudes. (This is indeed a remarkable result; the rank-2 problem entails a 
set of six parameters, but it has been shown that the set is reduced to that consisting 
of just three.) Therefore, the distribution function in both w and Q is given by 
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^-^^l^M.(u.,a)\'p, (C24) 



independent of the nucleon helicities. 

In another words, the spin density matrix has rank 1. The key ingredients for this result 
are that both spin-nonflip and spin-flip amplitudes harbor two resonances in D^ and P+ and 
that the production phase is mass-independent. It should be emphasized that the derivation 
given in this note does not depend on the existence of a good mass fit; it merely states that 
any fit with a mass-independent production phase is necessarily a rank-1 fit. Of course, the 
point is moot, if there exists no satisfactory fit in this model. 

rank condition derived 
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FIGURES 
FIG. 1. Experimental layout for E852. 

FIG. 2. Average acceptance vs. riTT~ effective mass. 

FIG. 3. Average acceptance vs. cos6 for different ryvr" effective mass regions. 

FIG. 4. Average acceptance vs. (p ^or different r/vr" effective mass regions. 

FIG. 5. Average acceptance vs. \t'\ integrated over all r/vr" effective masses. 

FIG. 6. Two-photon effective mass distribution. The central cross-hatched region shows the 
events which remain after SQUAW fitting. The shaded sidebands show the regions selected to 
estimate the background using Method 1 (see text). 

FIG. 7. Missing-massed squared distribution. The dashed histogram shows the distribution of 
events which remain after kinematic fitting. 

FIG. 8. Missing-mass squared vs. two-photon mass. The elliptical regions are used to estimate 
the background using Method 2 (see text). 

FIG. 9. Effective-mass distribution of the background estimated by Method 1. 

FIG. 10. Angular distribution of the background (Method 2) shown separately for (a.) the 
low-mass sideband and (b.) the high-mass sideband of the r]. 

FIG. 11. r]TT~ effective mass distribution. The shaded region is an estimate of the background 
using Method 2. 
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FIG. 12. Distribution of |t'| = \t — tmini- This experiment (solid dots) compared to a second 
experiment (open circles, see text) and to a Regge Pole fit. 

FIG. 13. Distributions of the acceptance-corrected cosine of the decay angle in the GJ frame 
for various effective mass selections. 

FIG. 14. Forward-backward asymmetry as a function of effective mass. The asymmetry 
= {F — B)/{F + B) where F(B) is the number of events for which the r/ decays in the forward 
(backward) hemisphere in the GJ frame. 

FIG. 15. Distributions of the Treiman-Yang angle (p in the GJ frame for various effective mass 
selections. 

FIG. 16. Effective mass distributions for: a.) the vr^p; and b.) the rjp systems. 

FIG. 17. Results of the partial wave amplitude analysis. Shown are a) the fitted intensity 
distributions for the D+ and b) the P+ partial waves, and c) their phase difference A$ . The 
range of values for the eight ambiguous solutions is shown by the central bar and the extent of the 
maximum error is shown by the error bars. Also shown as curves in a), b), and c) are the results of 
the mass dependent analysis described in the text. The lines in d) correspond to (1) the fitted !)+ 
Breit-Wigner phase, (2) the fitted P+ Breit-Wigner phase, (3) the fitted relative production phase 
0, and (4) the overall phase difference A<I>. 

FIG. 18. Results of the partial wave amplitude analysis. Shown are the fitted intensity distri- 
butions for the waves produced by unnatural-parity exchange. 

FIG. 19. Forward-backward asymmetry as a function of effective mass. Shown are: the total 
asymmetry in the data (closed circles); the predicted asymmetry from the PWA fit (open squares); 
the prediction of the fit for that part of the asymmetry due to natural-parity exchange (filled 
squares) ; and the prediction of the fit for that part of the asymmetry due to the unnatural-parity 
exchange waves (open circles). 
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FIG. 20. Experimental moments (open circles) shown with the predicted moments (open trian- 
gles) from the amplitude analysis. 

FIG. 21. Value of the log likelihood as a function of the number of P+ events in the PWA fit for 
all 8 ambiguous solutions. The inset shows a view with expanded scales. Because some solutions 
are very close to each other, not all 8 solutions are distinguishable on this figure. 

FIG. 22. Results of the partial wave amplitude analysis when the natural parity exchange m = 2 
amplitude is included. Shown are the fitted intensity distributions for a) the D+, b) the ?+ , and 
d) the D2+ partial waves. Shown in c) is the phase difference Ac/) between the D-(_ and the P+ 
partial waves. 

FIG. 23. Fitted P+ intensity and P+—D^ phase difference for the Monte Carlo sample generated 
with a pure D+ sample of a^(1320) events. 

FIG. 24. Comparison of the results of this amplitude analysis with the VES experiment. 

FIG. 25. Comparison of the results of this amplitude analysis with those of the KEK experiment. 

FIG. 26. The Tr+Tr^vr'^ effective mass distribution for events with the topology of three forward 
charged tracks, one recoil charged track, and two photon clusters consistent with a Tr*^. 

FIG. 27. The r/vr" effective mass distribution for the 77 -^ tt'^tt~tt^ event sample. 

FIG. 28. Comparison of the results of the amplitude analysis for the rj -^ tt^tt^tt^ (filled 
triangles) and the ij ^ 2'j (open circles) samples. 

FIG. 29. The fit results of the MDPWA (curves) and of one solution for the mass indepen- 
dent PWA (crosses) for the r]TT^ system: a) P+ , b) D-^- intensities and c) their relative phase 
A(p{P-^ — D+). Fig. la also shows the contributions of the 1 ^ signal intensity (1), the sum of the 
leakage and (signal - leakage) interference term (2) and the complete 1 ^ wave (3). 
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FIG. 30. Dependence of the 1 ^ signal parameters w^ (1), r^ (2) and the change in the 
-log(Likehhood) function relative to its minimum (3) on the leakage contribution TZi at (p = 80°. 
The black points are at the position of the likelihood extremum. 

FIG. 31. cos(a + Sa - 4) as a function of w from 1.2 to 1.6 GeV/c^ for a = 0°(o), a = 45°(+) 
andQ = 90°(n). 

FIG. 32. sin^ 5a (o), sin^ Jf, (+) and sin (5a sinJf, cos(a + Sa — S^) (□) as a function of w from 1.2 
to 1.6 GeV/c^, using a = 45°. 

FIG. 33. a^ sin^ 5a (o), b'^ sin^ 5b (+) and 2ab sin (5a sinSb cos(q + 5a — 5b) (n) as a function of 
w from 1.2 to 1.6 GeV/c^, where one has assumed that a = 1.0, b = 0.20 and a = 45°. 

FIG. 34. Angular distribution in cos^ as a function of w from 1.2 to 1.6 GeV/c^, where one 
has assumed that a = 1.0, b = 0.151 and a = 37.46°. 
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TABLES 



TABLE L Reduction in the data sample as a function of the cut type. 



Cut 


Remaining number of events 


Fraction removed (%) 


Number of triggers 


47 X 10^ 


- 


Topological and trigger cuts 


583,094 


98.8 


r] preselection (C.L. > 10""^) 


270,364 


53.6 


Removal of runs with LGD 


159,871 


40.9 


trigger processor failure 






LGD fiducial cut 


146,584 


8.3 


Photon-hadron distance cut 


145,710 


0.60 


Missing-mass-squared cut 


103,341 


29.1 


Confidence level cut 


85,888 


16.9 


A(/.<8° 


53,219 


38.0 


TPX2 cut 


49,113 


7.7 


Cut on position at EV/BV 


47,235 


3.8 


0.10 < \t\ < 0.95 (GeV/c)2 


42,676 


9.7 


0.98 < M{r]7r-) < 1.82 GeV/c^ 


38,272 


10.3 



TABLE IL Comparison of the results of E852 and the Crystal Barrel for the parameters of the 
J^'^ = 1 ^ resonance. 



Mass (MeV/c^ 



Width (MeV/c2 



E852 

Crystal Barrel 



1370 ±16 tlo 



1400 ±20 ±20 



385 ±40 



310 ±50 



+65 
-105 

+50 
-30 
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TABLE III. Comparison of the results of the PWA combined with a separate mass dependent 
fit (MDF) with those of the MDPWA with leakage. 



Meson 




Mass (MeV/c2) Width (MeV/c^) 


a^(1320) 


E852 (PWA+MDF) 
E852 (MDPWA) 


1317 ±1 ±2 127 ±2 ±2 
1313 ±1 119 ±2 


7rf(1400) 


E852 (PWA+MDF) 
E852 (MDPWA) 


1370 ±16 tlo 385 ±40 +^(^5 
1369 ±14 517 ±40 
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Figure 2 - Average acceptance vs. effective 
mass. 
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Figure 3 - Average acceptance vs. cos for 
different effective mass regions. 
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Figure 4 - Average acceptance vs. (]) for 
different effective mass regions. 
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Figure 6 - Two-photon effective mass 
distribution. The central cross-hatched region 
shows the events which remain after SQUAW 
fitting. The shaded sidebands show the 
regions selected to estimate the background 
using Method I (see text). 
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Figure 7 - Missing-massed squared 
distribution. The dashed histogram shows 
the distribution of events which remain 
after kinematic fitting. 
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Figure 8 - Missing-mass squared vs. 
two-photon mass. The elliptical regions 
are used to estimate the background using 
Method II (see text). 
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Figure 9 - Effective-mass distribution of the 
background estimated by Method I. 
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Figure 10 - Angular distribution of the 
background (Method I) shown separately for 
(a.) the low-mass sideband and (b.) the 

high-mass sideband of the r[. 
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Figure 1 1 - riTi: "effective mass distribution. The 
shaded region is an estimate of the background 
using Method 11. 



Figure 12 - Distribution of lt'l= It-t^^jj^ I. 
This experiment (solid dots) compared to a 
second experiment (open circles, see text) 
and to a Regge Pole fit. 
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Figure 13 - Distributions of the acceptance- 
corrected cosine of the decay angle in the GJ 
frame for various effective mass selections. 
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Figure 14 - Forward-backward asymmetry as a 
function of effective mass. The asymmetry 
=(F-B)/(F-i-B) where F(B) is the number of events 

for which the X] decays in the forward (backward) 
hemisphere in the GJ frame. 
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Figure 15 - Distributions of the Treiman- 

Yang angle (|) in the GJ frame for various 
effective mass selections. 



Figure 16 - Effective mass distributions 
for: a.) the 7l"p and b.) the Tjp systems. 
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Figure 17 - Results of the partial wave amplitude 
analysis. Shown are a) the fitted intensity 
distributions for the D^and b) the P^partial 

waves, and c) their phase difference A(|). The 
range of values for the eight ambiguous 
solutions is shown by the central bar and the 
extent of the maximum error is shown by the 
error bars. Also shown as curves in a), b), and 
c) are the results of the mass dependent analysis 
described in the text. The lines in d) 
correspond to (1) the fitted D+ Breit-Wigner 
phase, (2) the fitted P, Breit-Wigner phase, (3) 

the fitted relative production phase (|), and (4) the 

overall phase difference A4>. 
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Figure 1 8 - Results of the partial wave 
amplitude analysis. Shown are the fitted 
intensity distributions for the waves 
produced by unnatural-parity exchange. 
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Figure 19 - Forward-backward asymmetry as a 
function of effective mass. Shown are: the total 
asymmetry in the data (closed circles); the 
predicted asymmetry from the PWA fit (open 
squares); the prediction of the fit for that part of 
the asymmetry due to natural-parity exchange 
(filled squares); and the prediction of the fit for 
that part of the asymmetry due to the unnatural- 
parity exchange waves (open circles). 
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Figure 20 - Experimental moments (open circles) shown with the 
predicted moments (open triangles) from the amplitude analysis. 
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Figure 21 -Value of the log likelihood as a function 
of the number of P+ events in the PWA fit for all 8 
ambiguous solutions. The inset shows a view with 
expanded scales. Because some solutions are very 
close to each other, not all 8 solutions are 
distinguishable on this figure. 
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Figure 23 - Fitted P,. intensity and 1+ - D 
phase difference for the Monte Carlo 
sample generated with a pure D^. sample 
of a2(1320) events. 
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Figure 22 - Results of the partial wave amplitude 
analysis when the natural parity exchange m=2 
amplitude is included. Shown are the fitted 
intensity distributions for a) the D+, b) the I^, 
and d) the D2+ partial waves. Shown in c) is 

the phase difference A(|) between the E), and the 
P partial waves. 
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Figure 24 - Comparison of the results of this 
amplitude analysis with the VES experiment. 
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Figure 25 - Comparison of the results of this 
amplitude analysis with the KEK experiment. 









2500 






2000 


_ 




— 






o 


- 




^ 1500 


- 




CM 

O 


- 




1 1000 


- 


. 


> 

LU 






f\ 


500 


- 




^ \ 




, .1 


u 


v.^^ 



0.5 1.0 1.5 2.0 2.5 

M(7tVn°) (GeV/c^) 



3.0 



Figure 26 - The n n % effective mass distribution 
for events with the topology of three forward 
charged tracks, one recoil charged track, and two 

photon clusters consistent with a 7r. 
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Figure 27 - The r\ii effective mass distribution for 
the r\^ TZ n n event sample. 
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Figure 28 - Comparison of the results of the 
amplitude analysis for the r\^ n n n (filled 
triangles) and the T] ^ yy (open circles) samples. 
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Figure 29 - The fit results of the MDPWA (curves) 
and of one solution for the mass independent PWA 

(crosses) for the riJi system: a) P , b) D^. intensities 

and c) their relative phase A(|)(P - Q )• Fig. la also 
shows the contributions of the 1 ""*" signal intensity 
(1), the sum of the leakage and (signal - leakage) 
interference term (2) and the complete 1""*" wave (3). 
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Figure 30 - Dependence of the 1 signal 

parameters w j (curve 1), F i (curve 2) and 
the change in the -log(Likelihood) function 
relative to its minimum (curve 3) on the 

leakage contribution R at (|)"^ = 80°. The 
black points are at the position of the 
likelihood extremum. 
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Figure 31 - cos(a + 5^- d^) as a function of w 
from 1.2 to 1.6 GeV/c^ for a = 0° (diamonds), 
a = 45° (+) and a =90° (squares). 
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Figure 32 - sin 5 (diamonds), sin 5l(+) and 

sin 5j^sin 5, cos(a + 5 - 5, ) (squares) as a function 

of w from 1.2 to 1.6 GeV/c? using a = 45°. 
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Figure 33 - a siirSj(diamonds), b sin qX+), and 

2 a b sin S^sin S|^cos(a + 6^-613) (squares) as a 
function of w from 1.2 to 1.6 GeV/c^, where one 

has assumed that a = 1.0, b = 0.20 and a = 45°. 




Figure 34 - Angular distribution in cos(0) as a 
function of w from 1.2 to 1.6 GeV/c^, where one 

has assumed that a = 1.0, b = .151 and a = 37.46°. 



